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ABSTRACT 

As  a  first  step  towards  the  design  of  a  high  frequency,  high  force,  large  strain  shape  memory 
alloy(SMA)  actuator,  we  model  in  this  work  a  thermoelectrically  cooled  thin  SMA  layer  linear 
actuator.  The  SMA  is  subjected  to  cyclic  phase  transition  between  the  martensitic  and  austenitic 
phases  by  alternate  heating/cooling,  achieved  with  the  thermoelectric  Peltier  effect  of  a  pair  of 
P/N  semiconductors.  The  effect  of  variable  actuating  load  and  constant  load  applied  as  boundary 
conditions  on  the  SMA  actuator  are  considered.  The  thermomechanical  boundary  value  problem 
involves  strongly  coupled  thermal  and  mechanical  fields.  The  evolution  equations  for  the  field 
variables  are  integrated  using  the  fourth-order  Runge  Kutta  method  and  the  coupling  between  the 
fields  is  accounted  for  by  implementing  an  interative  scheme.  The  primary  parameters  of  interest 
in  this  work  are  the  frequency  response  and  evolution  of  the  variable  load.  The  performance  of 
the  actuator  is  compared  with  various  commercially  available  actuators  based  on  energy  conversion 
efficiencies  and  energy  output  per  unit  volume  of  active  material.  Results  of  the  analysis  indicate 
that  thin  SMA  layers(?^  6/i  thick)  under  partial  phase  transformation  are  capable  of  delivering 
frequencies  of  about  30  Hz  at  peak  stresses  of  about  145  MPa. 

1.  INTRODUCTION 

An  important  issue  related  to  the  dynamic  response  of  structures  is  the  suppression  of  unde¬ 
sirable  vibrations.  This  has  motivated  the  study  of  high  frequency,  low  force,  low  strain  actua¬ 
tors  using  piezoelectrics,  magnetostrictive  or  ferroelectric  mat erials(Takagi,  1990,  Wada  ei.aZ.,1990, 
Giurgiutiu  etai.,1995).  Of  more  recent  origin  is  the  study  of  certain  large  strain  actuators  which 
rely  on  the  imique  solid-solid  phase  transformation  of  shape  memory  alloys(SMA)(Wayman,1983). 
The  phase  transformation  is  highly  sensitive  to  an  applied  thermomechanical  loading  and  if  the 
imposed  stresses  are  reasonably  high,  significant  deformations  can  be  generated.  For  example,  the 
corresponding  actuation  forces  and  deformations  are  at  least  2  orders  of  magnitude  higher  than 
piezoelectrics.  The  shape  memory  effect  is  most  significant  in  Nitinol  which  is  composed  of  50 
at.wt.%  each  of  Nickel(Ni)  and  Titanium(Ti)(or  55%  and  45%  by  weight  of  Ni  and  Ti)  respectively. 
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The  phase  transformation  in  a  Ni-Ti  SMA  is  accompanied  by  a  significant  exchange  of  latent 
heat  with  the  environment.  The  rate  of  the  transformation  is  controlled  solely  by  the  time  rate  of 
heat  transfer  between  the  SMA  and  its  environment.  Conventional  heat  exchange  mechanisms  of 
resistive  heating  and  forced  convection  cooling  are  slow  and  hence  lead  to  a  low  frequency  response 
for  SMA  actuators(Bo  and  Lagoudas,1994).  A  novel  idea  to  increase  the  frequency  of  a  SMA 
actuator  was  put  forth  by  Lagoudas  and  Kinra(1993).  They  proposed  to  employ  the  heat  transfer 
due  to  the  thermoelectric  Peltier  effect(Domenciali,1954)  which  occurs  between  the  junction  of 
dissimilar  metals  due  to  the  flow  of  electric  current;  depending  on  the  current  direction,  one  junction 
becomes  cold  and  the  other  becomes  hot.  They  suggested  that  the  SMA  be  used  directly  as  one  of 
the  junctions,  which  can  then  be  heated  or  cooled  by  alternating  the  current  direction. 

To  demonstrate  this  principle,  consider  a  simple  SMA  actuator  shown  in  Fig.l.  A  thin  plate  of 
SMA  is  used  as  a  junction  by  sandwiching  it  between  a  pair  of  positively  doped(P)  and  negatively 
doped(N)  semiconductors.  The  other  junction  can  be  collectively  represented  by  heat  sinks(not 
shown)  positioned  at  the  ends  of  the  P  and  N.  If  the  SMA  layer  is  mechanically  constrained  at  the 
sides  and  the  phase  transformation  is  triggered  on  heating,  an  actuating  stress,  cr(t),  is  produced. 

Historically,  thermoelectric  coolers  have  been  used  primarily  as  refrigerators  and  thus  it  was 
their  steady  state  cooling  capability  that  was  of  interest.  Recently,  Thrasher  ef.aZ.(1992)  have 
theoretically  addressed  some  issues  related  to  the  study  of  the  transient  thermoelectric  phenomenon; 
a  comprehensive  study  however  appeared  to  be  lacking.  This  prompted  Bhattacharyya  et.aL{1995) 
to  extensively  study  the  transient  thermoelectric  problem  in  absence  of  actuation(no  mechanical 
constraints  on  the  SMA)  for  the  configuration  shown  in  Fig.l.  They  assumed  that  the  temperature 
field  is  one-dimensional  (along  the  z-axis),  and  found  that  due  to  the  high  thermal  diffusivity  of 
a  Ni-Ti  SMA,  the  temperature  gradients  in  a  thin  SMA  layer  are  negligible.  With  a  view  to 
investigating  periodic  solutions  of  thermal  cycling(alternate  heating  and  cooling),  Lagoudas  and 
Ding(1995)  reexamined  the  same  problem  by  assuming  at  the  outset  a  1-D  temperature  field  in 
the  P/N  elements  and  a  uniform  temperature  field  in  the  SMA.  In  this  paper,  we  extend  their 
study  to  the  issue  of  the  thermoelectromechanical  response  where  the  actuator  depicted  in  Fig.l 
is  subjected  to  a  mechanical  load  boundary  condition.  The  semiconductors  are  assumed  to  be 
stress-free  whereas  the  stress  state  in  the  SMA  is  taken  to  be  uniaxial  (along  the  x-direction). 

It  is  assumed  that  the  SMA  is  thermoelcLstic.  There  are  several  constitutive  models  which  have 
been  proposed  to  characterize  the  thermomechanical  response  of  a  polycrystalline  SMA  during 
phase  transformation.  A  non-exclusive  list  is  the  work  of  Tanaka(1986),  Patoor,  Eberhardt  and 
Berveiller(1987),  Liang  and  Rogers(1990,1991),  Sun  and  Hwang(1993a,b).  Recently,  Boyd  and 
Lagoudas(1995)  proposed  a  phenomenological  thermodynamic  theory  modeling  an  idealized  SMA 
response.  Lagoudas,  Bo  and  Bhattacharyya(1996)  extended  it  to  model  the  gradual  transformation 
of  SMA  polycrystals.  We  shall  use  the  latter  in  this  work.  The  resulting  actuation  stress  a{t)  of  the 
actuator  and  its  frequency  response  are  the  primary  parameters  of  interest.  The  frequency  range 
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studied  is  moderate  enough  so  that  inertial  effects  can  be  neglected.  Another  issue  of  interest  is  the 
coupling  between  the  thermal  and  mechanical  fields  in  the  SMA  during  the  phase  transformation. 

The  outline  of  this  paper  is  as  follows.  Sec.  2  presents  the  general  thermo  electromechanical 
governing  equations.  Sec, 3  presents  the  thermomechanical  boundary  value  problem  for  the  SMA 
layer  whereas  Sec.4  gives  the  thermoelectric  boundary  value  problem  for  the  P/N  semiconductors. 
Sec, 5  summarises  the  thermoelectromechanical  equations  pertaining  to  the  system  depicted  in  Fig.l, 
Sec.6  presents  the  results  of  a  parametric  study  and  is  followed  by  Sec. 7  on  the  conclusions. 

The  scientific  notation  used  in  the  text  will  be  now  described.  Vectors  will  be  denoted  with  a 
bold  face,  lower  case  Latin  letter(except  for  the  symbolic  vector  operator  del,  V)  ;  the  symbol  n 
will  be  reserved  for  a  unit  vector.  The  scalar  product  of  two  vectors,  x  and  n,  is  denoted  by  x.n 
and  it  has  the  evaluation  x,ni  in  terms  of  components  where  indicial  notation  has  been  used  and 
the  Einstein  summation  convention  is  implied.  On  the  other  hand,  their  tensor  product  is  denoted 
a:  (S)  n  ,  or  in  indicial  notation  as  XiUj,  resulting  in  a  second  order  tensor.  Magnitude  of  a  vector,x, 
will  be  denoted  cis  \x\  =  Second  order  tensors  will  be  denoted  with  a  bold  face,  lower 

case  Greek  letter  and  fourth  order  tensors  will  be  denoted  with  upper  case,  bold  face  Latin  letter. 
The  trace  of  a  second  order  tensor,  cr,  is  denoted  tr(cr)  and  in  indicial  notation  is  given  by  an. 
The  product  of  a  with  a  vector  v  is  av,  or  in  indicial  notation  aijvj.  The  product  of  two  second 
order  tensors,  a  and  A,  is  denoted  aX  and  its  evaluation  in  indicial  notation  is  aikXkj^  Their  scalar 
product  is  defined  as  a:X  —  tr(crA^) (superscript  T  indicates  transpose).  The  tensor  product  of  the 
two  second  order  tensors  is  denoted  o’  0  A,  or  in  components  as  aijXku  resulting  in  a  fourth  order 
tensor.  The  inverse  of  any  quantity,  will  be  denoted  as  ^4”^  and  the  material  derivative  will  be 
denoted  as  A.  While  the  indicial  notation  has  been  used  in  describing  the  mathematical  notation, 
the  various  tensor  operations  in  the  text  will  be  described  solely  by  symbolic  manipulation  of  the 
vector  and  tensor  quantities. 


2.  THE  THERMOELECTROMECHANICAL  GOVERNING  EQUATIONS 


2.1  The  Conservation  Laws 

Consider  a  given  mass  of  the  material  occupying  instantaneously  a  volume  V  surrounded  by 
the  surface  5.  We  shall  assume  that  body  forces  are  negligible  and  ignore  inertial  effects  (implying 
that  the  frequency  of  actuation  is  moderate;  we  shall  address  this  issue  in  Sec.6. 2  on  ^Tarametric 
Studies”).  For  a  particle  occupying  the  position  x  in  V  at  time  t,  the  equations  of  motion,  resulting 
from  the  conservation  of  linear  momentum,  reduce  in  the  absence  of  applied  body  forces  and  inertial 
effects,  to  the  equilibrium  equations 


div  =  0  , 


(2.1) 
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where  <t  =  cr{x^t)  is  Cauchy’s  stress  tensor.  The  angular  momentum  is  also  conserved,  leading 
to  the  symmetry  of  the  Cauchy  stress  tensor.  We  shall  restrict  the  analysis  to  the  study  of  SMA 
materials  for  which  its  mass  density  remains  unchanged  as  it  undergoes  a  phase  transformation 
between  its  two  crystalline  phases  of  austenite  and  maxtensite.  This  assumption  is  reasonable  as 
the  phase  transformation  in  SMAs  is  primarily  a  shear  process(Delaey  e^.aZ.,1974).  The  conservation 
of  mass  is  then  trivially  satisfied. 

The  conservation  of  energy  is  given  by 

PsmaU  =  -  V.g  +  r  +  cr:e,  (2.2) 

where  Psma  is  the  mass  density  of  the  SMA  material,  u  =  u{x^t)  is  the  specific  internal  energy 
rate(i.e.  internal  energy  rate  per  unit  mass),  q  =  g(x,  t)  is  the  energy  flux  vector(use  of  the  words 
“energy  flux”  instead  of  the  familiar  “heat  flux”  is  relevant  in  the  context  of  thermoelectricity,  to 
be  explained  later)  and  r  =  r(x,  t)  is  the  heat  source  per  unit  volume.  The  strain  rate  tensor  is 
denoted  €  s  €(x,  t),  given  for  small  strain  rates  and  rotations  by 

e  =  i[(V®tt)  +  (V®tif]  ,  (2.3) 

in  terms  of  the  rate  of  the  displacement  vector  it  =  n(x,t)(not  to  be  confused  with  the  scalar,  u, 
used  to  denote  the  specific  internal  energy).  This  is  a  reasonable  assumption  even  for  deformations 
during  phase  transformation  in  a  polycrystalline  SMA  where  the  total  strains  usually  do  not  exceed 
8%(Delaey  eta/., 1974). 

As  we  shall  see  later,  the  energy  flux  depends  on  the  current  density  vector,  j  =  j{x,t).  This 
implies  the  consideration  of  conservation  of  charge,  stated  as 

^  +  V.i  =  0.  (2.4) 

where  pc  =  pc{x,t)  is  the  net  charge  per  unit  volume.  While  electrical  transients  will  occur  as  a 
current  density  is  applied  at  t  =  0,  it  is  expected  that  the  time  interval  over  which  these  transients 
vanish  will  be  infinitesimal  compared  to  that  of  the  transients  in  the  thermal  and  mechanical  fields. 
Based  on  this  observation,  also  made  by  Lagoudcts  and  Ding(1995),  we  shall  assume  that  ^  =  0. 
We  then  have 


V.i  =  0.  (2.5)- 

2.2.  Thermoelectromechanical  Constitutive  Equations 

2.2.1  The  Shape  Memory  Alloy 

The  thermomechanical  constitutive  response  of  the  SMA  is  given  in  the  rate  form  by 

e  =  +  €«*  +  €«  =  Mo-  +  M&  +  c*'*  +  e‘  ,  (2.6) 
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where  €,  and  are  the  rates  of  the  total  strain,  elastic  strain,  thermal  strain  and  inelastic 

strain  due  to  phase  transformation,  respectively.  The  fourth  order  effective  elastic  compliance 
tensor  of  the  polycrystalline  SMA  material  is  denoted  as  M  and  will  be  defined  in  Eq.(2.9)  below. 
The  thermal  and  inelastic  strain  rates  are(Boyd  and  Lagoudas,1995) 

(2.7) 


where  is  the  coefficient  of  thermal  expansion  second  order  tensor,  T  is  the  temperature,  AT  = 
T  — To  with  To  being  a  reference  temperature  and  the  second  order  tensor  A  defines  the  “direction” 
in  which  the  inelastic  strains  develop  during  the  phase  transformation.  The  constituent  phases  of 
austenite(A)  and  martensite (M)  for  polycrystalline  SMA  are  assumed  to  have  isotropic  compliances 
and  coefficient  of  thermal  expansion  tensors.  Therefore 

Mi  =  ^I  +  ^iJ-S®S)  ,  a\>^=a\^S,  (2.8) 

Jbi  hi 

where  E**,  ui  and  are  the  Young’s  modulus,  Poisson’s  ratio  and  the  coefficient  of  thermal 
expansion,  respectively  for  the  ith  phase.  I  and  S  are  the  fourth  order  and  second  order  identity 
tensors  respectively. 

For  simplicity,  the  corresponding  effective  properties  of  the  SMA  are  approximated  by  the  rule 
of  mixtures  as(Boyd  and  Lagoudas,1995) 

M  =  M(0  =  Ma  +  ^ AM  and  a  =  =  cc^a  +  >  (2.9) 


where  AM  =  Mm  -  Ma  and  Aa^^  =  -  q^.  A  comparison  with  the  Mori-Tanaka  averag¬ 

ing  scheme,  as  discussed  by  Boyd  and  Lagoudas(1995)  for  a  polycrstalline  SMA  has  shown  that 
the  above  approximation  is  adequate.  The  tensor  A  is  assumed  to  have  the  following  representa- 
tion(Boyd  and  Lagoudas,1995) 

X  =  ,  ^>0  and  A  =  ,  ^<0,  (2.10) 

2  (Jg 

where  H  is  the  maximum  uniaxial  transformation  strain,  (jg  =  ^|cr' :  o-'j  ^  ,  <t'  =  o’  —  jtr (or)5  and 
is  the  inelastic  phase  transformation  strain  before  the  M  A  phase  transformation  commences, 
with  ^  •  The  evolution  of  the  volume  fraction  of  martensite  ^  is  given  by(Boyd 

and  Lagoudas,1995) 


f  Psma' 


dAuiO 


-1 


(2.11) 


where  V’  =  V’(0)/3  =  /^(^i CTyT).  The  explicit  forms  of  tp  and  f3  axe 


Psma  As 

7^^ 


AMcr  +  Aa'^AT  +  A 


(2.12) 
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In  the  above  expressions  As  =  sm  ~  sa  f'lid  Au(^)  =  UMiO  ~  the  average  specific 

entropy  and  u,(^)  is  the  internal  energy  of  the  ith  phase  at  the  reference  state  that  corresponds  to 
<r  =  0  and  T  =  Tq. 

Eqs.(2.7),(2.8)  and  (2.11)  in  Eq.(2.6)  result  in 

€  =  Mt&  +  Ot  r  ,  (2.13) 

where  Mt  is  the  tangent  compliance  fourth-order  tensor  and  is  the  tangent  coefficient  of 
thermal  expansion  second-order  tensor  during  the  transformation.These  are  given  by 

Mt  =  Mt (^,  tr,  T)  =  M(0  +  0-,  T)  ®  f3{^,  a,  T)  , 

oc^l^  =  a^^{^,<r,T)  =  a^HO  +  mm^,T)  ,  (2.14) 

In  absence  of  transformation(i.e.  ^  =  0),  we  set  P{^,ar,T)  =  0  aind  ip{Q  =  0  in  the  above 
expressions. 

During  phase  transformation,  Eq.(2.2)  becomes(Boyd  and  Lagoudas,1995) 

-  V.q  +  r  =  Cr -1- Ta*'*  :  d- -I- 7^  ,  (2.15) 

where  for  simplicity,  the  effective  heat  capacity,  C,  is  given  by  the  rule  of  mixtures(Boyd  and 
Lagoudas,1995)  as  C  =  C(0  =  Ca  +  ^(<^a/  —  Ca),  and  C,-  is  the  effective  heat  capacity  of  the  ith 
phase.  Note  that  in  absence  of  phase  transformation,i.e.  ^  =  0,  the  conservation  of  energy  equation 
(2.15)  reduces  to  the  well  known  form  for  a  thermoelastic  material(Boley  and  Weiner,1960).  It 
is  known  that  the  thermoelastic  coupling  term,  :  <r,  is  negligibly  small  during  quasistatic 

mechemical  loading(Boley  and  Weiner,1960).  With  phase  transformation  however,  additional  ther¬ 
momechanical  coupling  terms  are  introduced  through  the  term  7^(213  |  is  dependent  on  T  and  <7; 
see  Eq.(2.11)).  We  shall  investigate  this  coupling  in  Sec.5.4.  The  parameter,  7,  in  Eq.(2.15),  is 
given  by 

7(<r,  T)  =Ta  :  -b  psma^sT  -  11  ,  (2.16) 

where  11  is  the  thermodynamic  driving  force  for  the  phase  transformation  and  is  given  by 

n  =  :  AMcr  +  or  ;  (Aq'^AT  -b  A)  -b  psma^sT  -  Psma^u{0  ■  (2-17) 

Dming  the  phase  transformation,  11  attains  the  material  threshold  value  Y 

n=iy  for  ^<0,  y>o.  (2.18) 

The  method  of  determination  of  the  parameters  psma^s  ,  PsmaAu{()  in  Eq.(2.12)  and  Y  in  Eq.(2.18) 
from  experiments  axe  given  in  Sec.A  of  the  Appendix. 
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2.2.2  Thermoelectric  Materials 


There  is  no  phase  transformation  involved  in  a  thermoelectric  material.  Thus  Eq.(2.2)  simply 
becomes  the  heat  conduction  equation(Boley  and  Weiner, 1960) 

--V.g  +  r  =  Cr  ,  (2.19) 

for  a  stress-free  material(we  shall  return  to  this  point  in  Sec.4.1). 

In  order  to  give  the  constitutive  equations  of  thermoelectricity,  we  start  with  the  definition  of 
the  energy  flux  vector(Domenciali,1954) 

q  =  +  = -kVT  ,  =  (aT  + ,  (2.20) 

Qc 

where  the  vectors  and  are  the  thermal  and  electrical  components  of  the  energy  flux  vector 
respectively.  The  former  is  taken  to  relate  to  the  temperature  gradient  through  the  Fourier  law 
of  heat  conduction,  k  is  the  second  order  tensor  of  thermal  conductivity  and  a  is  the  second 
order  tensor  of  Seebeck  coefficients  (representing  the  phenomenon  of  thermoelectricity).  We  point 
out  that  the  temperature  T  used  in  the  statement  of  q^^  is  in  degrees  Kelvin,  as  required  by  the 
theory  of  thermoelectricity (Domenciali,  1954);  Qc  is  the  particle  charge  and  pL  is  the  electrochemical 
potential.  The  electrical  current  density,  jf,  is  (Domenciali, 1954) 

j  =  -p-H-V/i  +  a^VT)  ,  (2.21) 

Qc 

where  p  is  the  second  order  tensor  of  electrical  resistivity,  V/i  =  Vp  —  Qce  where  p  is  the  chemical 
potential  and  e  is  the  electric  field  vector.  Note  that  using  Eq.(2.20),  the  term  V.g  in  Eq.(2.19) 
becomes 

— V.g  =  K  :  (V  ®  VT)  +  div  k.VT  —  a.^'VT.j  —  (a  :  V  0  j)T  —  (div  ct^ .j)T 

_i(V/2.i-pV.j)  .(2.22) 

Qc 

Of  particular  interest  to  this  analysis  is  an  isotropic,  homogeneous  material  for  which  the  heat 
conduction  equation  is  simplified  considerably.  For  an  isotropic  material, 

p  z=z  pi  ^  K  =  Ki  ,  a  =  ai  .  (2,23) 

If  these  properties  are  assumed  to  be  spatially  uniform  and  insensitive  to  temperature  variations, 
then 

Vp  =  0  ,  V«  =  0  ,  Va  =  0.  (2.24) 

With  Eqs.(2.23)-(2.24),  div  k  =  Vk  =  0,  div  =  Va  =  0.  These  simplifications,  along  with 
Eqs.(2.5)  and  (2.21),  result  in 

- V.g  =  - V.g*'*  +  pIjP  =  /cV. VT  +  p|jf  .  (2.25) 
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It  is  seen  that  the  Joule  heating  term,  has  emerged  as  a  part  of  —  V.g.  The  parameter  r 

in  Eq.(2.15)  is  then  interpreted  as  collectively  representing  all  heat  sources  except  Joule  heating. 
In  the  next  section,  the  thermomechanical  boundary  value  problem(BVP)  for  the  SMA  layer  is 
presented. 


3-  THE  THERMOMECHANICAL  BOUNDARY  VALUE  PROBLEM  FOR  THE 

THIN  SMA  LAYER 

For  convenience,  we  present  the  mechanical  BVP  and  the  thermal  BVP  separately  in  Secs. 3.1 
and  3.2  respectively.  In  either  case,  as  the  corresponding  conservation  laws(i.e.governing  equations) 
have  been  already  introduced  in  Sec.2.1,  the  relevant  boundary  conditions  and  assumptions  to  sim¬ 
plify  the  problem  will  now  be  presented. 


3.1  The  Mechanical  Boundary  Value  Problem 


Based  on  the  specific  geometry  of  the  thin  SMA  layer  depicted  in  Fig.l,  we  assume  the  following 
boundary  conditions  based  on  a  cartesian  coordinate  system  with  its  origin  at  the  center  of  the 
SMA  layer. 


/+  w 
ern  dA  =  kg  (5o  —  n  at  x  =—  —  , 


(3.1) 


an  =  0  at 


(3.2) 


where  the  vector  n  is  taJken  as  the  outward  unit  normal  to  the  surface  of  the  SMA;  b,d  and  w  are 
the  dimensions  of  the  SMA,  as  indicated  in  Fig.l.  The  first  two  conditions  of  the  above  represent  a 
spring  loadled  boundary  condition  at  x  =—  ^5  ^^(t)  is  assumed  to  be  spatially  uniform,  the  spring 
has  a  stiffness  kg  and  an  initial  stretch  Jq-  The  remaining  surfaces  are  taken  to  be  traction-free  aind 
thus  includes  the  assumption  of  a  frictionless  interface  between  the  P/N.  The  temperature  and  the 
martensitic  volume  fraction  are  taken  as  initially  uniform  in  the  SMA  layer 


T  =  Tq  and  ^  =  ^0  t  =  0  . 


(3.3). 


We  shall  also  assume  that  the  initial  thermal  and  inelastic  strains  are 


=  0  and  c*  =  0  at  t  =  0  . 


(3.4) 


Denoting  Cxx  as  the  normal  component  of  a  in  the  x-direction,  it  is  seen  that  the  following  uniaxial 
stress  field 


axx  — ^^ia:(j/}  ^)  j 


(3.5) 
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with  all  other  components  of  or  being  zero,  satisfies  the  equilibrium  equation  (2.1)  and  the  traction- 
free  bovmdary  conditions  of  Eq.(3.2). 

At  this  point,  we  make  a  simplifying  assumption  that  the  temperature  field  is  spatially  uniform 
in  the  SMA  layer 

T  =  Tit)  .  (3.6) 

The  rationale  behind  the  above  assumption  will  be  discussed  in  Sec.3.2  on  “The  Thermal  Boundary 
Value  Problem”.  With  the  latter  of  Eq.(3.3),  and  Eqs.(3.5)-(3.6),  it  can  be  easily  inferred  from 
Eq.(2.11)  that  ^  =  ^(y,  2,t).  This  in  turn  implies  from  Eq.(2.14)  that  Mr  and  ar,  like  Cxx,  are 
independent  of  x  and  as  a  consequence,  we  conclude  from  Eq.(2.13)  that  e  =  eiy,z,t).  Therefore, 
in  absence  of  rigid  body  motion, 

iix  =  €xx  X  ,  (3.7) 


following  from  Eq.(2.3);  iix  is  the  rate  of  displacement  of  a  particle  in  the  x-direction  and  ixx  is 
the  normal  component  of  the  strain  rate,  €,  in  the  x-direction.  With  the  first  of  Eq.(3.1) (recall 
the  assumption  that  Us{t)  is  spatially  uniform),  we  note  that  Usit){—  ixx^  from  Eq.(3.7))  will  be 
spatially  uniform  only  if  ixx  is;  thus  by  Eq.(2.13),  <T(or  its  sole  nonvanishing  component  axx)  will 
necesscirily  have  to  be  spatially  uniform.  The  resulting  uniaxial  stress  field  is  denoted  as 


(Txx  =  cr(t)  .  (3.8) 

The  above  solution,  applied  to  the  second  of  Eq.(3.1),  results  in 

a(t)  =  ^(So-Us(t)).  (3.9) 

We  now  note  that  since  both  state  variables,  T(t)  and  cr(t),  are  spatially  uniform,  so  is  ^ (t)  (inferred 
from  Eq.(2.11)  and  the  second  of  Eq.(3.3)).  The  preceding  conclusions  also  imply  the  spatial 
uniformity  of  e  by  Eq,(2.13)  using  which,  with  (3.7)  and  the  rate  form  of  Eq.(3.9),  we  have 

a(t)  =  s(^,c7,T)f(t),  (3.10) 


where 


and 


=  - 


'2bdETi^,(7,T) 
.  w  ks 


ETi^,a,T)atlt{^,a,T)  , 


(3.11) 


ETi(,a,T)  =  [EiO-''+pH^,cr,T)]  '  and  a^^i^,<r,T)  =  . 


E{^)  is  the  effective  Young’s  modulus  of  the  SMA  layer,  and  P{^,a,T)  is  the  normal  component  of 
the  tensor  /3(^,  cr,T)  in  the  x-direction. 
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Thus  among  the  three  state  variables,  and  T,  the  evolution  of  the  first  two  are  given 
by  Eqs.(2.11)  and  (3.10)  respectively.  Evolution  of  the  temperature  field  will  follow  from  the 
conservation  of  energy  for  the  SMA  layer,  to  be  given  next  in  Sec. 3. 2. 

Before  we  go  on  to  the  next  section,  we  note  that  the  initial  values  of  the  state  variables  are 
needed  alongwith  their  evolution  equations  to  compute  their  total  values  at  a  time  t.  While  these 
have  been  prescribed  for  T  and  ^  by  Eq.(3.3),  the  initial  value  of  stress,  cr(0),  will  follow  from  the 
boimdary  conditions  of  Eqs.(3.1),(3.2)  and  the  solution  given  by  (3.8)  respectively.  At  t  =  0,  using 
Eq.(2.6)  for  purely  elastic  behavior  and  with  Eq.(3.4),  we  obtain 

where  the  spatially  uniform  normal  component,  €xx,  of  e  in  the  x-direction  is  denoted  as  e{t). 
Combining  the  above  with  the  total  form  of  Eq.(3.7)  and  the  second  of  Eq.(3.1),  we  have 


(7(0)  = 


bd  + 


2E{0) 


(3.13) 


3.2  The  Thermal  Boundary  Value  Problem 

We  now  recall  Eq.(3.6),  wherein  we  made  the  assumption  of  an  uniform  temperature  field  in 
the  SMA  thin  layer  (thickness  of  the  order  of  mm).  This  assumption  has  its  origin  in  the  transient 
thermoelectric  problem  addressed  by  Bhattacharyya  eLa/.(1995)  in  absence  of  actuation  forces 
where  a  1-D  variation(along  the  z-axis)  of  temperature  in  the  P/SMA/N  geometry  of  Fig.l  was 
assumed.  The  assumption  of  a  1-D  temperature  field  was  verified  by  the  finite  element  method 
and  was  found  to  be  an  excellent  one.  The  method  of  separation  of  variables  was  invoked  to  solve 
the  1-D  problem  analytically.  It  was  found  that,  as  a  consequence  of  the  high  thermal  diffusivity 
of  a  Ni-Ti  SMA,  the  temperature  field  was  virtually  uniform  in  the  SMA  layer.  Lagoudas  and 
Ding(1995),  in  an  effort  to  find  periodic  solutions  of  a  cyclic  thermal  field,  addressed  the  same 
problem  but  assumed  an  uniform  temperature  field  in  the  SMA  at  the  outset.  The  analytical 
solution  resulted  in  an  integro-differential  equation  for  the  SMA  temperature  which  they  solved  by  a 
finite  difference  scheme.  Prompted  by  these  considerations,  we  also  assume  an  uniform  temperature  - 
field  in  the  SMA;  our  formulation  then  is  an  extension  of  Lagoudas  and  Ding’s(1995)  work  to  the 
thermoelectromechanical  case. 

We  assume  that  the  electric  current  density  vector  is  directed  along  the  z-direction,  such  that 

j  =  J{t)n,  .  (3.14) 

The  portion  of  the  surface  of  the  SMA  in  interfacial  contact  with  the  P/N  elements  is  defined 
as  Sj;  the  remainder  of  the  surface  is  denoted  as  Se*  The  thermal  boundary  condition  on  5/  is 
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defined  as 


=  q  on  Sj  :  z  =—  ^  ,  (3.15) 

where  the  heat  flux  q  will  come  from  the  interface  condition  between  the  SMA  and  P/N  and  will 
be  made  explicit  in  Sec. 4.  A  convective  boundary  condition  is  specified  for  the  remaining  surfaces 
of  the  SMA.  Therefore 

q.n  = -h(T{t) -To)  on  Se  x  =- ^  ,y  ^  ,  (3.16) 


where  h  is  the  convection  coefficient  and  To  is  the  ambient  temperature.  Note  that  the  initial 
condition  for  the  temperature  is  given  by  the  first  of  Eq.(3.3). 

We  now  write  Eq.(2.15)  alongwith  Eq.(2.25)  in  an  integral  form  and  use  Eqs.(3.6),  (3.10),  (3.16) 
and  (3.14),  to  get 


f  qdS-2h(l  +  -){T{t)-To)+pJ{t)^=\c{0  +  7i<r,T) 
Vsma  Js, 


I  Psma 


^Au(0 


-1 


m+ 


f  Psma~ 


5Au(^) 


-1 


s{^,a,T)}f{t)  ,(3.17) 


where  the  volume  of  the  SMA  layer  is  denoted  as  Vsma  and  the  joule  heating  is  accounted  for  by 
setting  r  =  pJ{t)^  in  Eq.(2.15),  where  p  is  the  resistivity  of  the  SMA.  We  now  have  three  evolution 
equations  for  the  state  variables,  a  and  T  given  by  Eqs.(2.11),  (3.10)  and  (3.17)  and  these  can 
be  solved  once  the  heat  flux  q  in  Eq.(3.17)  is  specified.  This  is  done  in  Sec.4  next.  We  point  out 
here  that  all  quantities  pertaining  to  the  SMA  have  been  denoted  without  any  subscripts  (except  its 
density  psma,  volume,  Vsma  and  mass,  Msma)-  In  the  next  section  where  the  thermoelectric  problem 
for  the  P/N  semiconductors  will  be  discussed,  all  quantities  pertaining  to  the  P/N  semiconductors 
will  be  labeled  with  the  subscript  P  or  N  respectively. 


4.  THE  THERMOELECTRIC  BOUNDARY  VALUE  PROBLEM  FOR  THE 

SEMICONDUCTORS 


As  in  the  previous  section  for  the  SMA  layer,  we  first  present  the  boundary  conditions  for.^ 
the  semiconductors  and  certain  relevant  assumptions.  Any  simplification  in  the  conservation  laws, 
already  presented  in  Sec.2.2,  will  then  be  discussed.  We  shall  assume  that  the  tractions  are  con¬ 
tinuous  at  the  interfaces  of  the  SMA  with  P/N.  Thus,  by  Eq.(3.2),  we  conclude  that  the  interface 
is  traction-free.  Assuming  that  the  remaining  surfaces  of  the  P/N  elements  are  also  traction-free, 
the  mechanical  boundary  conditions  for  P/N  become 

<Tn  =  0  on  Si  {i  =  P,N),  (4.1) 
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where  5,-  is  the  surface  bounding  the  ith  phase.  The  solution  for  the  stress  field  in  the  P/N  elements 
then  follows  simply  as 


o-  =  0  .  (4.2) 

The  temperature  gradient  and  the  current  density  vector  are  taken  to  be  unidirectional  (along 
the  z-axis) (Bhattacharyya  ef.a/.(1995);  see  discussion  at  the  beginning  of  Sec. 3. 2).  These  conditions 
are  summarized  as 


T  =  Ti{z,t)  and  j  =  J{t)n,  (i  =  P,N)  .  (4.3) 

An  isothermal  boundary  condition  is  assumed  at  the  end  surfaces 

{i  =  P,N).  (4.4) 

The  interfaces  between  the  SMA  with  P  and  N  are  taken  to  be  thermally  perfect,  so  that  the 
temperature  and  the  energy  fiux  are  continuous.  These  are  stated  respectively  as 

T{t)  =  Ti  and  q.n  =  q,.n  at  z  =- ^  (i  =  P,N)  ,  (4.5) 

where  n  is  an  outward  unit  normal  on  the  surface  of  P/N  elements.  We  shall  use  the  latter  boundary 
condition  shortly  to  give  a  discussion  of  the  Peltier  effect.  A  convective  boundai'y  condition  is 
assumed  for  the  remaining  surfaces 

g,-.n  =  -h{Ti{z, t)  -Tq)  on  x  ^  ,y  ^  {i  -  P, N)  .  (4.6) 

The  initial  condition  for  the  temperature  is  taken  as 

T.(z,0)=To  ii  =  P,N).  (4.7) 

The  conservation  of  energy  equation  will  follow  from  Eq.(2.19)  and  Eq.(2.25)  as 

KiV.VTi  +  PiJ^ it)  +  r  =  Citi  (i  =  P,N)  .  (4.8) 

The  convection  boundary  condition  given  by  Eq.(4.6)  for  the  ith  element  is  included  approximately^ 
as  a  source  term(refer  Bhattacharyya  et.al. ,1995);  thus  r  =  — 2/i  {Ti{z,  t)  — To).  In  addition, 

using  Eq.(4.3),  Eq.(4.8)  then  becomes 

/I  IN  PfP‘ 

i^i-^{z,t)  +  piJ{tf-2hy--\-—^{Ti{z,t)-TQ)=Ci-^{z,t),  {i  =  P,N),  t>0,  (4.9) 

where  ki,  pi  and  Ci  are  the  thermal  conductivity,  the  electrical  resistivity  and  the  heat  capacity  of 
the  ith  element  respectively. 
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We  now  give  a  discussion  of  the  Peltier  effect  based  on  the  interface  condition  for  the  energy 
flux,  given  by  the  latter  of  Eq.(4.5).  Starting  with  the  SMA/P  interface  and  assuming,  i.e.  p,  = 
/2p(Bhattacharyya  ei.aZ., 1995),  the  second  of  Eq.(4.5)  with  Eq.(2.20)  is  written  for  the  SMA/P 
interface  as 

(9*^  -  =  {Qp  -  9®^)  (9*^  “  9p^)  =  (“p  “  a)T{t)j.n;  ,  (4.10) 

where  n  =  -n^.  The  Seebeck  coefficients  of  the  P  and  SMA  are  ctp  and  a  respectively.  The  jump 
in  the  heat  flux,  denoted  by  Eq.(4,10),  represents  the  thermoelectric  Peltier  effect.  Note  that  this 
effect  is  driven  by  the  electric  current  and  influences  the  temperature  gradients.  This  is  exactly 
opposite  to  the  thermoelectric  Seebeck  effect,  based  on  which  thermocouples  operate.  In  this  latter 
ca^e,  if  the  junctions  between  two  dissimilar  metals  forming  a  closed  circuit  are  maintained  at 
two  different  temperatures,  an  electric  current  will  flow.  Both  effects  are  controlled  by  the  same 
parameter,  i.e.  the  Seebeck  coefficient(Harman  and  Honig,1967). 

Returning  to  Eq.(4.10),  the  jump  in  the  heat  flux  indicates  that  the  Peltier  effect  has  manifested 
itself  as  a  heat  source/sink  at  the  interface.  For  commercially  available  Bismuth  Telluride(Bi- 
Te)  semiconductors,  ap  =  2.15  x  10"''^yAr"’^(Melcor,1992)  whereas  a  =  1.2  x  lO^^FAT^^  (Jackson 
ei.aZ.,1972)  for  a  Ni-Ti  SMA  and  ap-a  >  0.  We  note  that  T{t)  has  to  be  in  absolute  units(Domenciali,1954), 
i.e.r(t)  >  0  for  t  >  0.  Examining  the  right-hand  side  of  Eq.(4.10),  it  is  clear  that  when  the  current 
density  vector  is  directed  from  the  N  to  the  P(j.n^  >  0),  we  have  >  0*  This  implies 

that  the  heat  lost  by  the  SMA  at  the  SMA/P  interface  is  in  excess  of  that  which  is  conducted  into 
the  P.  The  difference,  for  the  chosen  current  direction,  represents  the  Peltier  effect  acting  as  a  heat 
sink  at  the  SMA/P  interface.  Should  the  current  density  be  reversed,  it  is  easy  to  deduce  that  the 
Peltier  effect  now  acts  as  a  heat  source. 

Let  us  now  consider  the  situation  at  the  SMA/N  interface.  With  Eq.(2.20),  the  latter  of  Eq.(4.5) 
is  written  for  this  interface  as 

(9"*  -  Qn)  =  {Qu  -  9'')  ^  (9*''  -  9n'‘)  =  («„  -  a)T{t)j.n,  ,  (4.11) 

where  n  =  and  assuming  p  =  For  Bi-Te  semiconductors  =  — ap(Melcor,1992);  thus 
an  —  QJ  <  0.  From  the  above  equation,  we  see  that  if  the  current  is  directed  from  the  N  to  the 
>  0),  then  .n^  <  0.  Thus,  at  the  SMA/N  interface,  heat  lost  by  the  SMA  is  in  - 

excess  of  that  which  is  conducted  into  the  N.  Therefore,  the  Peltier  effect  at  this  interface  has  a 
similar  effect  as  at  the  SMA/P  interface  for  the  same  current  direction,i.e.  it  acts  as  a  heat  sink. 

When  the  current  is  reversed(P  to  N),  the  effect  is  opposite,  but  again  similar  at  both  interfaces, 
i.e.  the  Peltier  effect  acts  as  a  heat  source  at  both  interfaces. 

We  note  that  when  the  current  is  directed  from  the  N  to  the  P,  the  heat  lost  to  the  Peltier  heat 
sink  at  both  interfaces  has  the  potential  to  cool  the  SMA  inspite  of  the  fact  that  electrical  current 
is  being  passed  through  the  SMA  (providing  a  Joule  heat  source).  It  was  shown  by  Bhattacharyya 
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ei.aZ.(  1995)  in  their  theoretical  predictions  and  comparisons  with  experiments,  that  for  a  moderate 
current  density  and  with  a  given  geometry  of  P/SMA/N,  the  SMA  can  be  indeed  cooled  by  passing 
the  ciirrent  from  the  N  to  the  R  The  theoretical  predictions  however  show  that  there  is  a  critical 
current  density  above  which  it  is  not  possible  to  cool  the  material (Lagoudas  and  Ding, 1995)  due 
to  a  substantial  increase  of  Joule  heating.  If,  however,  the  current  was  directed  from  P  to  iV",  the 
Peltier  heat  source  reinforces  the  Joule  heating;  thus,  for  this  current  direction,  heating  occurs  at 
all  current  densities.  While  semiconductors  of  various  metals/alloys  are  available,  commercially 
available  BhTe  has  one  of  the  highest  magnitudes  of  ap  and  an  around  room  temperature  and  are 
thus  capable  of  delivering  the  maximum  thermoelectric  effect  (maximum  heating/cooling)  (Harman 
and  Honig,1967,Pollock,1991) 

Bhattachaxyya  et,aL{1995)  and  Lagoudas  and  Ding(1995)  invoked  symmetry  considerations 
based  on  certain  assumptions  to  simplify  the  mathematical  problem.  We  shall  do  the  same.  The 
assumptions  are  listed  below. 

1.  All  material  properties  of  the  P  and  N  Bi-Te  semiconductors  are  assumed  identical(Melcor,1992) 
except  for  their  Seebeck  coefficients  which  have  equal  magnitudes  but  opposite  signs.  Thus 

Kp  —  Kn  j  Pp  ^  pn  5  ^p  “  ^  ^p  ~  *  (^•^^) 

2.  Since  for  Bi-Te  semiconductors  and  Ni-Ti  SMAs,  ap  =  —an  ^  18a,  we  assume 

ap  =  —an  >>  a  — ap  —  a  ^  ap  ^  an'-  a  ^  an. 

The  above  assumption  has  been  tested  in  a  finite  element  computation  of  the  thermoelectric 
problem.  With  the  values  of  the  Seebeck  coefficients  for  P/SMA/N,  the  evolution  of  the  tempera¬ 
ture  in  the  SMA  was  computed  and  compared  with  the  case  when  a  for  the  SMA  was  set  to  zero. 
The  difference  in  the  computed  temperature  over  a  time  evolution  of  30  sec  was  less  than  0.5%. 
For  the  Seebeck  coefficients  of  Bi-Te  and  Ni-Ti,  the  assumption  is  thus  a  sound  one. 

With  the  above  assumptions,  the  1-D  temperature  distribution  is  symmetrical  about  the  plane, 

2  =  0.  This  was  shown  analytically  by  Lagoudas  and  Ding(1995).  Thus 

Tp{z,t)-Tn{-z,t)  ,  ^<z<^  +  dp  ,  t>0.  (4.13) 

With  Eq.(4.13),  only  half  of  the  P/SMA/N  assemblage  needs  to  be  considered;  we  therefore  shall 
focus  on  the  P  element  only,  i.e.  in  the  remainder  of  the  text,  Eq.(4.9)  will  be  considered  only  for 
the  P  element.  We  shall  now  combine  the  equations  developed  in  Secs.3  and  4  and  present  them 
in  their  final  form  in  the  next  section. 

5.  SOLUTION  OF  THE  COUPLED  THERMOELECTROMECHANICAL 
PROBLEM  AND  THE  ENERGY  CONVERSION  EFFICIENCY  OF  THE 

ACTUATOR 
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5.1  The  Solution  of  the  Coupled  Thermomechanical  Problem 

Recall  that  the  heat  flux  quantity  q  in  the  first  term  of  Eq.(3.17)  needs  to  be  specified  from  the 
interface  conditions  of  the  SMA  with  P/N.  With  the  symmetry  in  the  temperature  field  given  by 
Eq.(4.13),  it  is  sufficient  to  derive  the  heat  flux  ^(defined  first  in  Eq.(3.15))  either  from  Eq.(4.10) 
or  (4.11)  respectively,  with  the  help  of  the  second  of  Eq.(2.20)  and  Eq.(4.3).  It  is 

q  =  •  (5-^) 

Eq.(3.17)  reduces  to 

dAu(^)' 

psma 


[a‘'*(e)T(0  +  7(^,T) 

The  evolution  equations  (2.11)  and  (3.10)  for  ^  and  a  can  be  solved  alongwith  the  above  equation 

dT  fi 

once  the  heat  flux  term,  is  known.  This  implies  that  a  fourth  equation  is  needed;  it 

is  Eq.  (4.9)  (Due  to  the  symmetry  in  the  temperature  field,  we  shall  only  consider  the  equation  for 
the  P  element). 

In  absence  of  a  stress  field(when  only  Eqs.(4.9)  and  (5.2)  will  remain,  with  s(^,c7,T)  =  0  in 
the  latter),  the  thermoelectromechanical  problem  reduces  to  the  one  addressed  by  Lagoudas  and 
Ding(1995).  They  show  that  it  is  possible,  due  to  the  linearity  of  the  thermal  problem  in  the  P 
material,  to  derive  an  integro-differential  equation  solely  in  terms  of  the  temperature  of  the  SMA; 
the  derived  equation  does  not  involve  the  temperature  field  of  the  P/N  elements.  They  solved  the 
equation  by  the  finite  difference  scheme. 

For  the  current  problem,  the  methodology  of  Lagoudas  and  Ding(1995)  is  used  to  derive  the 
following  differential  equation  from  Eqs.(4.9)  and  (5.2) 

^{t)  =  S{t)  -  uiT{t)  .  t>0,  (5.3) 

where  the  temperature  of  the  SMA,  T(t),  can  be  solved  from  the  above  equation(and  with  the 
initial  condition,  given  by  the  first  of  Eq.(3.3))  once  S{t)  is  found  from  the  solution  of  the  integro- 
differential  equation 

f  G{t  -  r)SiT)dT  +  m(^ ,  a,  T)S{t)  +  [u2{t)  -  a,  T)ui]  f  6""^ (‘-’■)5(r)dT 
Jo  Jo 

=  F{t)  -  To  [U2{t)  -  ^(^,a,T)i.i]e-"“  .  (5.4) 

The  parameters  G{t),  i'2it),  and  F{t)  are  given  in  the  Appendix.  Note  that  the 

coupling  parameter  j{a,T)  enters  into  Eq.(5.4)  through  the  expression  for  fj.{^,a',T){see  the  third 
of  Eq.(A.ll)  in  the  Appendix). 


dAuiO 


—  X 


-  2/1  (i  +  i)  (T(t)  -  To)  +  pJ{tf  =  I  C{0  + 


(5.2) 
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To  summarise,  Eqs.(2.11),  (3.10)  and  (5.3)  correspond  to  the  evolution  equation  for  the  marten¬ 
sitic  volume  fraction,  the  conservation  of  linear  momentum  and  the  conservation  of  energy  respec¬ 
tively.  This  system  of  first-order  nonlinear  ordinary  differential  equations  alongwith  the  initial 
conditions  given  by  Eqs.(3.3)  and  (3.13)  and  supplemented  by  Eq.(5.4)  need  to  be  solved  for  a 
ajid  r.  The  computational  procedure  is  given  in  Sec  6.1. 


5.2  Energy  Conversion  Efficiency  of  the  Actuator 

Before  we  move  into  the  numerical  scheme,  it  is  desirable  to  characterize  the  performance  of  the 
thermoelectric  SMA  actuator.  Normally,  two  different  measures  are  used(Giurgiutiu  ei.ai.,1995). 
These  are: 

(a)  The  energy  output  per  unit  volume  of  the  active  material  which  is  stored  in  the  actu¬ 
ated  structure(the  spring  in  our  case)  during  the  transformation  of  the  SMA  from  martensite  into 
austenite  during  heating.  If  the  transformation  starts  at  t  the  energy  stored  during  the 

transformation  is  given  by 

In  case  of  a  very  compliant  spring(i.e. implying  a  constant  load  boundary  condition)  with  some  ini¬ 
tial  stress,  Eq.(5.5)  is  the  work  done  per  unit  volume  of  the  active  material  in  displacing  a  constant 
load  during  the  M  A  transformation  by  heating. 


(b)  The  associated  energy  conversion  efficiency  in  converting  the  electrical  energy  input  to  the 
actuator  into  the  energy  stored  in  the  spring  during  the  M  — >  A  transformation.  This  is  defined  as 


r){t)  = 


(5.6) 


where  we  define  W^{t)  as  the  electrical  energy  input  per  unit  volume  of  the  P/SMA/N  assembly. 
It  is  given  as 

=  Ive.jdVdT,  (5.7) 

where  V  is  the  volume  of  the  P/SMA/N  assembly.  Using  Eq.(2.21),  the  assumption  V/i  =  0  in  the 
P/SMA/N  and  the  last  of  Eq.(2.23),  Eq.(5.7)  becomes 


where 


W%t) 


W^rna{t)Vsma  +  2W^{t)Vp 

Vsma  +  2Vp 


W!^,{t)d+2W^{t)dp 

d  "H  2dp 


(5.8) 


Wtma{t)  =  pf^j\r)dr  and  W;(<)  =  ^  j\To -T{r))J{T)dT  +  pp  j‘^{T)di 


(5.9) 
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At  this  point,  we  note  that  a  thermoelectric  device  has  been  historically  used  eis  a  refrigerator 
and  to  characterize  its  performance,  the  definition  of  “Coefficient  of  Performance”  (COP)  (Nag, 1981) 
has  been  used.  We  shall  now  provide  a  point  of  contact  between  the  notion  of  actuator  efficiency 
and  the  COP  of  a  thermoelectric  device. 

Whenever  a  thermoelectric  device  has  been  used  as  a  refrigerator,  its  steady  state  “cooling” 
capability  has  been  of  interest  and  the  COP  has  been  defined  in  that  context.  In  our  case,  it  is 
the  transient  heat  transfer  that  we  are  interested  in;  moreover  the  thermoelectric  device  is  not  only 
being  used  as  a  refrigerator  but  also  as  a  heat  pump.  As  we  compute  the  efficiency  of  the  actuator 
when  the  device  is  being  used  as  a  heat  pump,  we  shall  define  the  COP  in  the  context  of  transient 
heating  only.  The  COP  is  defined  as 


COP{t)  = 


^sma[Q{^)  Q(^o)] 

Weit)  -  Weito) 


(5.10) 


where  Msma  is  the  mass  of  the  SMA  layer  and  0(t)  is  the  specific  heat  input(heat  input  per  imit 
mass)  into  the  SMA(defined  in  Eq.(5.11)).  The  numerator  in  the  above  equation  represents  the 
total  thermal  energy  brought  into  the  SMA  layer  at  a  temperature  r(t)(with  the  heat  sinks  at  To) 
during  the  transformation.  Q{i)  is 


Q{t)  =  [  Q{T)dT  , 

Jo 


(5.11) 


where  Q(t)  follows  from  the  right  hand  side  of  Eq.(3.17)  as 

,  -1 


W  =  — ^C^(0  +  7(c^,T) 


I  Psn 


dAtx(^) 

'~W 


^>(0+ 


a^H0nt)  +  j{<r,T) 


/  Psma’ 


8Au{^) 


ae 


-1 


si^,cr,T)\f{t).  (5.12) 


The  COP  of  a  corresponding  Carnot  heat  pump(i.e.  the  maximum  attainable  COP)  is 


COPcarnotit)  = 


(5.13) 


With  Eqs.(5.10)  and  (5.13),  we  can  write  Eq.(5.6)  as 


rj{t)  = 


Msma[Qit)  -  Q(to)] 


COPjt) 

COPcarnotit 


X  COPcarnotit  • 


(5.14) 


Recognizing  ’M,ma[Q(i)--Q(to)]  ~  ^  efficiency  by  which  the  SMA  layer  converts  heat  into 

work  and  co^frnlit)  ~  as  the  COP  of  the  thermoelectric  device  relative  to  a  Caxnot 

device,  we  then  have 


T]it  —  ^sma(t)  X  COPrdit^  X  COP carnoti^')  •  (5.15) 
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In  addition,  it  is  reasonable  to  identify  r]sma(i)  x  COPcamotit)  =  Vcamotii)  as  the  efficiency  of  the 
SMA  layer  when  operated  with  a  Carnot  heat  pump.  Thus 

vii)  =  Vcarnot{t)  X  COPrel{t)  .  (5.16) 

Therefore,  rjcamotii)  is  the  maximum  attainable  efficiency  of  the  actuator  whereas  the  efficiency, 
r]{t),  of  the  “real”  actuator  is  controlled  by  the  COPrelif)  of  th®  thermoelectric  device. 

5.3  Resume  of  the  Nondimensionalized  Solution 


While  the  theoretical  development  was  based  on  dimensional  quantities,  it  is  preferable  to 
present  the  numerical  results  in  a  non-dimensional  form.  We  thus  obtain  the  nondimensional  coun¬ 
terparts  of  physical  quantities  by  the  application  of  the  Buckingham  tt  theorem(Buckingham,1914), 
which  are  listed  below.  Note  that  for  the  thermoelectric  problem,  the  relevant  nondimensional 
groups  were  originally  introduced  by  Bhattacharyya  et.a/.(1995). 


Normalized  absolute  temperature:  T  =  — —  1  , 

To 

a  .  _  (CpdlY^ 

Normalized  stress:  a  =  —  ,  Normalized  time:  i  =  t  ^  , 

Em  \  Kp  J  ^ 

Normalized  convection  coefficient:  /i  =  2h  f  ~  ^  —  ,  Normalized  current:  J  =  , 

\b  w J  Kp  y  f^pTo 


Normalized  Seebeck  coefficient:  dip  —  api 


I  To 

Kppp 


,  Normalized  coefficient  of  thermal  expansion:  a  =  aTo  , 


Normalized  stress  influence  coefficient:  D  = 

k. 


D  -  E 

— — J—  ,  Thermomechanical  Coupling  Coefficient:  Etc  =  'rr  ri  ’ 

Em  f  2o  Cp 


Normalized  Spring  Stiffness:  = 


Em'w 


,  Normalized  initial  stress:  a'(O)  = 


o'(Q) 
Em  ’ 


Relative  Heat  capacity  of  SMA:  C  =  ,  Relative  resistivity  of  SMA:  p  =  —  , 

Cp  pp 

E 

Relative  Young’s  modulus  of  SMA:  E  =  ^7—  , 

Em 

Relative  dimensions  of  the  SMA  and  the  P/N:  —  ,6=  —  ,cL=  —  . 

w  w  ^  w 

Using  the  nondimensional  groups  listed  in  Eq.(5.17),  we  present  a  resume  of  the  final  set  of 
equations  and  the  associated  initial  conditions,  all  in  a  nondimensionalised  form.  Eq.(3.10)  for  the 
conservation  of  linear  momentum  in  the  SMA  reduces  to 


(5.17) 


where 


5(e,a,f)  =  - 


=  s{^,a,T)T{i)  , 
J^T(^,d-,r)ar(^,a,r) 


2^ma  I  2 

ks 


(5.18) 


(5.19) 
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where  Eri^.a^f)  =  ET{^,a,T)/EM,  ari^.a^f)  =  Toari^.a^T)  and  the  normalized  volume  of 
the  SMA  is  Vsma  =  Eq.(2.11)  corresponding  to  the  evolution  equation  for  the  volume  fraction 

of  martensite  becomes 


where 

^(0  =  ^0 


f  Psma  ■ 


5Au(^) 


V-d)  and  P{^,a,T)-EM 


(5.20) 


f  psma 


dAu{^) 


-1 


/3(^,<T,r),(5.21) 


where  the  time  derivatives  in  Eq.(5.20)  are  with  respect  to  the  nondimensional  time,  t.  Finally, 
Eq.(5.3)  for  the  conservation  of  energy  in  the  SMA  reduces  to 


^(t)  =  S{t)  -  Pi [1  +  f  (t)]  ,  i  >  0  . 


(5.22) 


where  5(t)  =  and  Pi  is  given  in  Eqs.(A.12)  of  the  Appendix.  The  solution  to  the  above 

differential  equation  caji  be  found  once  the  parameter  S{t)  is  solved  from  the  nondimensionalised 
version  of  Eq.(5.4),  given  as 

f  G(t-f)5(f)df +  /Z(^,d,f)5(t)  +  [P2(i)  - /i(^,o-,r)Pi]  /  e“‘'‘^‘“’'^5(f)df 

Jo  Jo 

=  F{i)  -  [p2(t)  -  .  (5.23) 

The  parameters  G(t),  ^2(^5  and  F{t)  are  given  in  Eqs.(A.12)  of  the  Appendix.  The 

initial  conditions  for  the  nondimensional  stress,  the  martensitic  volume  fraction  and  the  nondimen¬ 
sional  temperature,  follow  from  Eqs.(3.3)  and  (3.13),  as 


m 


^0 


Vjfpia  I  1 

A:,  ^  2£(0) 


,  e(o)  =  6  ,  f(o)  =  o. 


(5.24) 


The  nondimensional  energy  stored  during  the  transformation  follows  from  Eq.(5.5)  as 

The  energy  efficiency(see  Eq.(5.16))  is  intrinsically  nondimensional,  therefore 

^(t)  =  7?(t)  . 

Eqs.(5.18)-(5.26)  constitute  the  nondimensionalised  solution. 


(5.25) 


(5.26) 


5.4  Thermomechanical  Coupling 
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Before  the  issue  of  thermomechanical  coupling  during  the  phase  transformation  is  addressed, 
we  list  some  typical  parameters  for  a  Ni-Ti  SMA.  These  are(Jackson  etal.,1972) 

=  23®C  ,  mJ  =  1°C  ,  4  =  29°C  ,  A}  =  51°C  , 

Ha  =  0.148  J/mm^  ,  H  =  8%  ,  Sp  =  3.5  , 

Ea  =  70x  10^  MPa  ,  Em  =  30  x  10^  MPa  ,  va  =  =  1/3  , 

qa  =  11  X  10-V°<^  >  “M  =  6.6  X  10“V°<^  .  (5-27) 

where  and  are  the  stress-free  martensitic  start  and  finish  temperatmes,  A®  and  A®  are  the 
austenite  start  and  finish  temperatures,  Ha  is  the  magnitude  of  the  latent  heat,  H  is  the  maximum 
uniaxial  inelastic  transformation  strain  and  Sp  is  the  standard  deviation  of  the  normal  distribution 
p{T)  defined  in  Eq.(A.6)  of  the  Appendix.  Ei  and  ui  are  the  Young’s  modulus  and  the  Poisson’s 
ratio  of  the  ith  phase(i=A,M).  With  these  values,  we  compute 

D  =  6.187  MPaj^C  and  Y  =  6.93  MPa  ,  (5.28) 

from  Eq.(A.lO)  in  the  Appendix.  The  thermoelectric  properties  of  the  SMA  and  the  P  element 
are  summarized  in  the  table  below.  Note  that  the  thermoelectric  properties  of  the  N  element  will 
follow  from  Eq.(4.12). 


a{V/K) 

/^(J/(mm.s-K) 

p(Q-mm) 

C  (J/(mm5  —  K)  ) 

SMA 

1.2  X  10-5 

2.2  X  10-2 

6.3242  X  10-“* 

2.12  xlO-5 

P 

2.15  X  lO"'* 

1.63  X  10-5 

1.15  X  10-2 

4.35  xlO-5 

Table  1.  Material  Constants  of  SMA  and  Semiconductor  elements. 

The  geometric  parameters  are  assumed  to  be 

b  —  w  =  4  mm  ,  dp  =  2  mm  ,  d  =  0.5  mm  .  (5.29) 

The  convection  coefficient  is  taken  as 

h  =  1  X  l0~^J/{mm^sK)  ^  .  (5.30) 

The  above  value  of  h  was  found  suitable  by  Bhattacharyya  et.a/.(1995)  to  model  free  convection 
conditions.  The  initial  values  of  the  stress,  martensitic  volume  fraction  and  the  temperature  are 
taken  as 

a(0)  =  0  ,  ^(0)  =  1  ,  r(0)  =  297°  K  .  (5.31) 

The  spring  stiffness  and  the  current  density  are  chosen  as 

kg  =  3  KN/mm  ,  J  =  3  Amps/mm^  .  (5.32) 
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Based  on  the  above,  we  now  give  the  nondimensional  parameters.  Note  that  among  those  listed 
in  Eq.(5.17),  the  evolution  of  the  nondimensional  temperature  and  the  stress,  T  and  d,  with  the 
nondimensional  time,  t,  are  of  interest.  The  paremeters  d.  Etc  3-nd  E  depend  on  the  state  variable, 

^  and  hence  will  follow  from  the  solution.  The  remaining  nondimensional  parameters  therefore  are 

h  =  2.454  X  10"^  ,  Op  =  0.856  ,  D  =  6.125  x  lO'^  ,  C  =  0.487  , 
p  =  0.055  ,  5  =  1  ,  dp  =  0.5  ,  (5.33) 

whereas 

J  =  0.925  ,  d(0)  =  0  ,  ^5  =  0.025  ,  d  =  0,125  .  (5.34) 

Among  the  nondimensional  parameters,  only  those  in  Eq.(5.34)  will  be  treated  as  variable  param¬ 
eters  in  the  numerical  studies;  hence,  their  numerical  values  will  be  mentioned  in  the  inset  of  the 
figures.  The  computational  scheme  based  on  which  the  numerical  results  were  generated  has  been 
discussed  in  detail  in  Sec. 6.1. 

In  order  to  understand  the  thermomechanical  coupling  in  the  heat  conduction  equation  Eqs.(5.22)- 
(5.23),  it  is  perhaps  more  useful  to  start  from  the  nondimensional  form  of  Eq. (2.15)  (from  which 
Eqs.(5.22)-(5.23))  were  derived).  It  is 

-V.g  +  f  =  +(l  +  f)a‘'‘:^  +  7^  <  (5.35) 

Etc 

where  among  the  nondimensional  parameters  involved  in  the  above  equation,  those  not  listed  in 
C  C  cP 

Eq.(5.17)  are  V.q  =  ^  and 7  =  'y/E!\^{\  note  that  the  rate  terms  in  Eq.(5.35)  are 

with  respect  to  the  nondimensional  time.  Our  intention  here  is  to  uncover  the  effect  of  the  coupling 
term  7^.  Towards  that  end,  we  look  at  the  evolution  of  T  vs.  i  and  a  vs.  i  given  in  Fig.2a  and 
Fig.2b  respectively  for  two  cases:  (a)  the  line  1  corresponds  to  omitting  the  term,  7^  from  Eq.(5.35) 
and  (b)  the  line  2  corresponds  to  retaining  the  term  7^.  In  both  cases,  the  SMA,  starting  from 
an  initially  fully  martensitic  condition,  is  heated  until  it  is  fully  converted  into  austenite(^  ==  0), 
following  which  the  current  direction  is  reversed  to  cool  the  SMA  layer  until  the  austenite  fully 
converts  to  martensite.  The  effect  of  the  phase  transformation  is  seen  to  slow  the  heating  and  the 
cooling  process. 

With  the  material  parameters  given  in  the  beginning  of  this  section,  it  can  be  shown  that 
7  =  jIEm  =  (TcrAa-f  pAsoT  —  II) /Em  ^  -'Ha/ Em  where  Ha  is  the  magnitude  of  the  latent  heat 
of  transformation(i7a  >  0).  Thus,  the  term  7^(or  ^^0  represents  the  latent  heat  exchanged 
during  a  phase  transformation.  When  ^  <  0 (during  heating),  this  term  acts  as  a  heat  sink  and 
indicates  the  endothermic  nature  of  the  M  A  transformation;  the  effect  is  to  slow  the  temperature 
increase,  evident  from  a  comparison  of  the  heating  portion  of  curve  2  with  that  of  curve  1  in  Fig.2a. 
During  the  A  M  transformation  upon  cooling,  the  term  acts  as  a  heat  source  and  is 

indicative  of  the  exothermic  nature  of  the  A  M  phase  transformation.  This  results  in  a  slower 
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cooling  process;  compare  cooling  portions  of  curves  2  and  1.  These  trends  in  the  evolution  rate  are 
also  reflected  in  the  evolution  of  a,  in  Fig. 2b. 

We  also  examine  the  issue  by  estimating  the  relative  magnitudes  of  the  terms  on  the  right  hand 
side  of  Eq. (5.35).  As  the  SMA  is  heated  from  ^  =  1  to  0(refer  the  line  2  in  Figs. 2a  and  b),  the  change 
in  T  and  a  are  approximately  AT  =  0.22  and  Aa  =  7.5  x  10"”^  and  the  change  in  ^  is  A^  =  —1.  Over 
the  time  range  that  this  transformation  has  taken  place,  the  cumulative  change  in  the  terms  on  the 
right  hand  side  of  Eq.(5.35),  will  be  of  the  order  of  ^^AT,  (1  +  T)a^^Aa  and  7A^,  respectively; 
these  are  estimated  to  be  1.979  x  10”^,  2.99  x  10”^  and  4.93  x  10"^,  respectively.  It  is  seen  then 
that  the  thermoeleistic  coupling  term(the  second  term  in  the  list)  is  relatively  very  small  whereas 
the  heat  capacity  term(the  first  term)  and  the  coupling  term  due  to  phase  transformation(the  last 
term)  are  of  the  same  order  of  magnitude,  with  the  latter  being  the  more  predominant  one.  This 
suggests  that  while  the  thermoelastic  coupling  term  may  still  be  omitted,  the  one  due  to  the  phase 
transformation,  7^,  must  be  retained.  This  has  not  always  been  the  practice  in  the  literature 
where  both  terms  have  been  omitted  to  enable  decoupling  of  the  heat  conduction  equation  from 
the  mechanical  problem  for  numerical  simplicity(Brinson, 1994, Oberaigner, 1995). 

Before  closing,  we  further  assess  the  relative  importance  of  the  stress  rate  term  and  the  tem¬ 
perature  rate  term  constituting  ^  in  Eq.(5.20).  The  ratio  of  the  net  evolution  of  each  term  is  of  the 
order  of  — £'A/Ad-/(roDAf );  this  ratio  follows  from  Eqs.(5.21),  (91b)  and  the  discussion  in  Sec. A 
of  the  Appendix.  Its  numerical  value  turns  out  to  be  around  —0.56.  Therefore,  the  temperature 
rate  term  and  the  stress  rate  term  are  of  the  same  order.  The  combined  effect  of  both  terms  is 
responsible  for  the  thermomechanical  coupling  during  phase  transformation. 


6.  THE  COMPUTATIONAL  SCHEME  AND  PARAMETRIC  STUDIES 
6.1  The  Computational  Scheme 

The  system  of  evolution  equations  (5.18),  (5.20)  and  (5.22)  will  be  solved  for  d,  ^  and  T 
incrementally  in  time,  using  the  fourth  order  Runge  Kutta  method  in  a  nonlinear  manner.  However, 
since  the  coefficients  of  these  equations  themselves  depend  on  d,  ^  and  T,  an  iterative  procedure 
additionally  is  required  at  every  time  step.  Thermomechanically  coupled  problems  have  been 
addressed  in  the  past,  see  Allen(1991)  and  Armero  and  Simo(1991).  For  our  specific  problem,  we 
outline  the  numerical  approach  adopted.  We  shall  first  discuss  the  implementation  of  the  Runge- 
Kutta  method  and  then  the  iterative  scheme. 

If  At  is  the  duration  of  a  time  increment,  we  define 

Cn  =  ^(nAt),dn  =  d(nAt)  and  r„  =  r(nAt)  ,  (6.1) 

for  each  positive  integer  n;  given  the  values  and  fn-i  at  the  (n-l)th  time  step,  we  need 

to  compute  the  values  and  Tn  at  the  nth  time  step.  This  is  done  by  first  writing  Eq.(5.18) 
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as 


an  -a„_l  =  5(r,a*,r)(f„  -Tn-l)  , 


and  Eq.(5.20)  as 

Cn  -  ^n-1  =  ^(r)(fn  -  fn-l)  +  ;S(r, a*, T*)(an  -  d„_i)  , 
where  initially,  we  assign 

C  =  ^n-i  ,  or*  =  an-i  and  T*  =  r„_i  . 


(6.2) 


(6.3) 


(6.4) 


The  temperature  increment  is  found  by  integrating  Eq.(5.22)  by  the  fourth-order  Runge  Kutta 
method(Numerical  Recipes,  1986);  this  approach  is  somewhat  more  accurate  than  the  finite  differ¬ 
ence  scheme  used  by  Lagoudas  and  Ding(1995).  Thus 


^  ^  ki  k2  kz  k^ 

Tn-T„_l  =  -  +  y  +  y  +  y  , 


(6.5) 


where 


S(£)  is  given  by 


ki  ^  At[S{tn^l)  —  l/l  —  UiT{tn^l)]  , 
k2  =  At 


At.  ^  _  ki 

S{tn-i  +  — )  -ui-  iyiT{tn-i)  - 


At^  ^  ko 

S{tn-l  +  — )  -Ui-  UiT{tn^i) 


A:3  =  At 

k4  =  At  [5(tn)  Ulf{in^i)  -  Uik^]  , 


Sit)  =  Sd 


^  [G(t  -  f)  +  [i>2(t)  -  /2(r,^*,T*)t>i]  df 


(6.6) 


SD  =  Fn-  [92{l)-pi{C,r,T)v,]e-^Fi 

«“1  rJfcAt 

-y^Skl  [G{i-  f)  +  [92i-piC,&\f-)9i]]  e-"‘‘df  . 

k=l 


and 

-j  rtiAi 

=  F(nAi)  ,  Sn  =  5n(nAt)  ,Gn  =  ^  G{i)  dt ,  Rn  =  5(nAt)  .  (6.7) 

J{n-l)£i.i 

Note  that  we  solve  Eq.(6.5)  for  fn  and  thereby  Eqs.(6.2),  followed  by  (6.3).  We  now  have  the  first 
set  of  trial  values  for  a„,  and  Tn.  If  the  difference  |r„  —  r*|  >  5(where  5  is  a  small  number),  then 
we  set  ^*,  a*  and  T*  equal  to  the  newly  computed  updated  values  and  then  recalculate  a  second  set 


23 


of  trial  values  from  Eqs.(6.5),  (6.2)  and  (6.3).  This  procedure  is  repeated  until  |r„  —  f*]  <  5;  the 
convergence  for  an  and  Tn  is  then  achieved.  We  mention  at  this  point  that  while  the  convergence 
criterion  has  been  based  solely  on  fn,  our  numerical  computations  confirm  excellent  simultaneous 
convergence  for  ^  and  an- 

6.2  Parametric  Studies 

We  recall  that  in  Sec. 5.4  on  “The  Thermomechanical  Coupling”,  it  was  mentioned  that  a, 
Etc  2md  E  depend  on  the  state  variables  ^  and  will  follow  from  the  solution  of  the  problem; 
the  parameters  are  given  in  Eq.(5.33)  and  will  be  assumed  to  remain  unchanged  in  the  rest  of 
the  parametric  study.  Thus,  the  influence  of  the  parameters,  J,  a(0),  ks  and  d  on  T  and  a  as  a 
function  of  t  is  of  interest.  While  the  parametric  study  will  be  primarily  based  on  non-dimensional 
parameters,  we  shall  end  with  a  sample  case  based  on  dimensional  parameters  in  order  to  get  a 
feel  for  the  physical  quantities  involved.  We  note  at  this  point  that  the  frequency  of  cycling  is 
dependent  on  the  amount  of  austenite  that  is  allowed  to  form(or  the  amount  of  martensite  that 
remains)  at  the  end  of  the  heating  cycle.  We  shall  thus  introduce  the  parameter  ^res  to  denote  the 
residual  martensite  volume  fraction  at  the  end  of  the  heating  cycle  and  assume  ^res  =  0.8  in  the 
following  analysis. 

Initially,  we  shall  investigate  the  influence  of  the  parameters  ks  and  d’(O);  the  former  corresponds 
to  the  stiffness  of  the  actuator  controlled  structure  and  the  latter  is  the  initial  stress  exerted  on  the 
actuator  due  to  an  initial  structural  deformation.  We  shall  fix  J  =  0.925  and  d  =  0.1 25  (the  values 
initially  assumed  in  Eq.(5.34))  and  the  evolution  of  T  and  a  shall  be  studied  for  a  range  of  kg  and 
a(0). 

In  order  to  simulate  a  constant  load  and  a  variable  load  boundary  condition,  we  consider  two  dif¬ 
ferent  spring  stiffnesses,  kg  =  8.3x  10~^  and  2.5  x  10“^.  Fig.3a  gives  the  nondimensional  temperature 
for  both  stiffnesses  with  the  same  initial  nondimensional  stress,  a-(O)  =  1.67  x  10^^(corresponding 
to  50  MPa  with  the  dimensional  parameters  given  in  Sec. 5. 4).  It  is  seen  from  Fig. 3b  that  the 
nondimensional  stress(and  by  implication,  the  dimensional  one)  does  not  change  during  the  course 
of  the  transformation  for  the  more  compliant  spring;  the  value  of  kg  =  8.3  x  10”^  is  thus  sufficient 
to  simulate  a  constant  load  boundary  condition. 

We  now  return  to  Fig.3a.  It  is  seen  that  the  two  curves  for  T  at  the  two  different  spring 
stiffnesses  are  virtually  indistinguishable  up  to  a  certain  point,  following  which  these  separate.  The 
onset  of  the  M  -4  A  transformation  is  almost  identical  in  both  cases,  A  slight  difference  occurs 
because  while  both  start  out  with  an  identical  stress  state,  the  thermal  expansion  of  the  SMA  on 
heating(and  prior  to  the  onset  of  the  M  -4  A  transformation)  alters  the  stress  state  differently  in  the 
two  cases.  However,  since  the  thermal  strains  axe  usually  very  small,  the  stress  states  in  both  ca^es 
are  only  marginally  different.  Recalling  that  the  transformation  temperatures  are  highly  sensitive  to 
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the  magnitude  of  the  uniaxial  stress  (increasing  with  stress  and  vice-versa),  “almost”  identical  stress 
states  imply  “almost”  identical  transformation  temperatures.  This  implies  that  the  onset  is  almost 
identical  in  both  cases.  However,  with  the  contraction  due  to  phase  traiisformation(substantial  in 
comparison  to  the  thermal  expansion),  the  subsequent  evolution  of  the  stress  in  the  two  cases  is 
substantially  different,  as  is  obvious  from  Fig. 3b.  The  increasing  stress  state  in  the  stiffer  spring 
raises  the  transformation  finish  temperature  and  thus  delays  the  completion  of  the  transformation, 
in  comparison  to  the  compliant  spring;  hence  with  a  stiffer  spring,  the  frequency  of  cycling  will 
always  be  lower  than  the  one  with  a  more  compliant  spring(with  the  same  initial  stress  in  both 
cases) . 

The  nondimensional  energy  output  and  the  energy  conversion  efficiency  rj  defined  in 

Eqs.(5.25)  and  (5.26)  are  computed  at  the  end  of  the  M  A  transformation(corresponding  to 
Crcs  =  0.8).  Their  values  at  various  initial  stress  levels  and  spring  stiffnesses  have  been  shown  in 
Figs.4a  and  4b(in  these  figures,  the  computed  data  points  have  been  indicated  with  open  circles 
and  then  connected  by  straight  lines;  such  a  procedure  hcis  been  also  adopted  in  Figs. 4c  through 
7).  It  is  seen  in  Fig. 4a  that  at  a  given  initial  stress  level,  an  increasing  spring  stiffness  leads  to  a 
greater  work  output.  This  can  be  anticipated.  With  a  higher  spring  stiffness,  the  stress  rate  will 
also  be  higher.  Thus  at  the  end  of  the  transformation,  while  identical  inelastic  strains  axe  recovered 
in  all  cases,  the  stress  corresponding  to  a  stiffer  spring  will  be  higher  and  thus  will  lead  to  a  greater 
work  output.  However,  the  energy  conversion  efficiency(see  Fig.4b)  decreases  with  a  stiffer  spring 
at  a  given  stress  level;  this  physically  happens  because  with  a  higher  stress  rate,  the  material  needs 
to  go  to  higher  temperatures  to  complete  the  transformation  leading  to  the  enhanced  consumption 
of  electrical  energy.  For  comparison,  we  also  give  the  corresponding  values  of  rjcamoti  defined  prior 
to  Eq.(5.16).  It  is  seen  that  the  maximum  achievable  efficiency  is  an  order  of  magnitude  higher 
than  the  real  efficiency  77.  This  implies(see  Eq.(5.16))  that  the  COP  of  a  thermoelectric  heat  pump 
is  an  order  of  magnitude  less  efficient  compared  to  a  Carnot  thermoelectric  heat  pump. 

At  this  point,  it  is  desirable  to  compare  these  performance  measures  to  other  commercially 
available  actuators.  Giurgiutiu  e^.aZ.(1995)  have  reported  a  comparative  study  of  ferroelectric, 
electrostrictive  and  magnetostrictive  actuators.  The  energy  output  per  unit  volume  is  the  highest 
for  an  electrostrictive  actuator,  =  11.913  J/m^^  whereas  a  magnetostrictive  actuator  has  the 
highest  reported  energy  conversion  efficiency  of  77  =  67.1  %.  For  the  thermoelectric  SMA  actua¬ 
tor,  we  find  that  =  2.884  x  10“^,  implying  that  =  86.52  x  10^  J/m^,  and 

7]  =  1.35%,  attainable  in  a  test  case  corresponding  to  the  parameters  used  to  generate  the  full 
transformation (^res  =  0)  cyclic  response  curve(line  2)  in  Fig.2.  Therefore,  the  SMA  thermoelec¬ 
tric  actuator,  in  comparison  with  the  actuators  studied  by  Giurgiutiu  e^.aZ.(1995),  has  a  very  high 
energy  output  per  unit  volume  of  active  material(four  orders  of  magnitude)  and  a  very  low  en¬ 
ergy  conversion  efficiency.  These  results  can  be  anticipated  since  SMAs  have  very  high  recoverable 
strains  (assumed  8%  in  our  case)  compared  to  about  0.143%  for  the  electrostrictive  material(from 
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Table  3  of  Giurgiutiu  ei.a/.,1995);  thus  the  specific  energy  output  for  the  SMA  actuator  is  ex¬ 
pected  to  be  considerably  higher.  On  the  other  hand,  a  SMA  relies  on  the  conversion  of  heat  into 
mechanical  energy  whereas  electrostrictive  and  magnetostrictive  materials  convert  electrical  and 
magnetic  energies  directly  into  mechanical  energy.  Thermal  energy  being  of  a  lower  grade  com¬ 
pared  to  electrical  and  magnetic  energies,  the  conversion  of  the  former  into  useful  mechanical  work 
is  accomplished  with  a  lower  efficiency  as  compared  to  the  conversion  of  the  latter. 

While  it  was  instructive  to  compare  the  energy  output  and  the  energy  conversion  efficiency  of  the 

SMA  actuator  with  that  of  other  commercially  available  actuators,  the  key  issue  for  SMA  actuators 

is  to  address  the  question  of  actuator  frequency.  In  the  remaining  part  of  the  nondimensional 

parametric  study,  we  shall  confine  ourselves  to  the  issue  of  frequency  response.  We  therefore 

turn  to  the  nondimensional  frequency  of  actuation,  depicted  in  Fig.5.  If  a;  is  the  frequency  of 

C  (P 

actuation(number  of  cycles  per  second),  then  Q  =  It  is  seen  that  at  both  initial  stress  levels, 

the  frequency  decreases  with  increasing  spring  stifFness(this  was  first  noted  when  discussing  Fig.Sa). 
In  order  to  get  an  idea  about  the  dimensional  frequency,  we  choose  a  sample  case;  i.e.  ks  ==  0.025 
and  cr(0)  =  3.33  x  10~^,  for  which  u  =  13.9  from  Fig.5.  The  corresponding  dimensional  frequency 
with  the  parameters  given  in  Sec. 5. 4  is  then  1.3  Hz. 

In  Fig.6,  we  fix  a(0)  =  3.33  x  10“^  and  the  ratio  d/ks  =  5.  Bound  by  this  condition,  we  study 
the  influence  on  Q  due  to  a  change  in  d  and  k^.  This  is  done  at  two  different  levels  of  the  current 
density,  J.  It  is  seen  that  a  decreasing  thickness  d(and  a  consistent  change  in  at  a  given  value  of 
the  current  J,  leads  to  a  higher  nondimensional  frequency,  a).  Now  let  us  examine  the  dimensional 
frequency,  a;,  for  a  test  case:  J  =  1.23  and  d  =  2.5  x  10“^  for  which  Q  =  102.6,  from  Fig.6.  With  the 
relevant  parameters  given  in  Sec. 5.4,  these  nondimensional  quantities  translate  to  a;  =  9.616  Hz  at 
J  =  4  Anvpsjmrr?  ^  d  =  10  fi.  It  can  thus  be  concluded  that  very  thin  SMA  actuators  have  relatively 
fast  response.  As  an  aside,  we  note  that  the  two  curves  intersect  each  other  around  d  >  0.11.  This 
is  because  for  thicker  SMAs  during  cooling,  the  Joule  heating  begins  to  dominate  over  the  Peltier 
effect  when  the  current  density  exceeds  a  certain  critical  value (Lagoudas  and  Ding,1995);thus  the 
cooling  rate  decreases  or  the  time  needed  to  cool  down  to  a  certain  temperature  increases  and 
reduces  the  frequency  of  actuation.  When  the  SMA  thickness  is  reduced  below  d  «  0.11,  this 
critical  value  exceeds  the  normalized  current  densities  shown  in  the  figure.  Thus  below  d  =  0.11,  (D 
increases  as  J  increases  from  0.925  to  1.23. 

Until  now,  we  have  not  discussed  the  influence  of  material  parameters  on  the  frequency.  To 
be  specific,  it  is  well  known  that  the  hysteresis  in  the  stress-free  transformation  temperatures  in  a 
NiTi  SMA  are  highly  sensitive  to  alloying  elements  (Jackson  ei.a/.,1972)  and  thus,  it  is  a  parameter 
of  interest.  The  influence  of  the  hysteresis  on  the  frequency  is  now  studied. 

Recall  that  the  hysteresis  in  the  stress-free  transformation  temperatures  is  responsible  for  the 
dissipation  during  a  phase  transformation  in  the  SMA.  In  the  constitutive  theory  for  SMAs  we 
have  used  here,  a  reference  to  the  work  of  Boyd  and  Lagoudas(1995)  shows  that  the  term,  /Hd^, 
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represents  the  cumulative  dissipation  during  the  phase  transformation  and  over  a  complete  cycle 
of  the  M  A  and  A  M  transformations,  it  is  straightforward  to  show  from  Eq.(2.18)  that  the 
net  dissipation  is  2y ;  y  is  thus  a  measure  of  dissipation.  We  shall  now  establish  the  link  between 
the  parameter  y  and  the  transformation  temperature  hysteresis.  For  simplicity,  we  assume  that 
the  spread  between  the  start  and  finish  temperatures  during  both  transformations  is  identical, i.e. 

AJ  -  -M]  -4  AI-MJ  =  a)-  M°  .  (6.8) 

We  now  define  the  trajisformation  temperature  hysteresis  as  the  parameter  ATam  =  -^5  ~  = 

A®  —  Mg,  such  that  A®  =  +  ATam  and  A°  =  Mg  +  ATam-  With  these  definitions,  we  rewrite 

y  from  the  second  of  Eq.(A.lO)  in  the  Appendix  as 

1  r  fMj  rMUATAM  1 

y  =  j-HD  /  Tp{T)dT  -  /  Tp{T)dT  .  (6.9) 

Notice  that  if  ATam  =  0,  then  y  =  0;  this  implies  that  in  absence  of  hysteresis,  the  dissipation 
vanishes.  The  above  equation  then  links  the  measure  of  dissipation  to  the  hysteresis  in  the  trans¬ 
formation  temperatures.  With  respect  to  the  temperatures  in  Eq.(5.27),  the  first  of  Eq.(6.8)  shows 
that  A® —  A3  =  M3  —  =  22^C  whereas  ATam  =  28^C  from  Eq.(5.27).  For  the  parametric  study, 

we  keep  the  former  fixed  and  change  the  hysteresis  ATam  in  the  range  20  <  ATam  ^  28  (reducing 
ATam  below  20  raises  issues  regarding  minor  loops(Mayergoyz,1991),  and  is  beyond  the  scope  of 
this  paper).  It  is  seen  from  Fig. 7  that  a  decrease  in  the  hysteresis  results  in  a  substantial  increase 
in  frequency.  With  ATam  =  20®C,  Q  =  143;  the  corresponding  dimensional  frequency  based  on  the 
dimensional  parameters  in  Sec.5.4  is  a;  =  13.4.  If  however  we  set  J  =  6  Amps/mm^,  then  we  have 
dp  =  1.33  mm,  b  ^  w  —  2.67  mm,  d  =  6.67/x  and  the  dimensional  frequency  becomes  a;  =  30.15  Hz; 
the  evolution  of  the  dimensional  temperature  and  stress  are  shown  in  Fig.8. 

We  now  return  to  Eq.(2.1)  where  the  inertial  effects  were  ignored  and  the  conservation  of  linear 
momentum  reduced  to  the  equilibrium  equation.  We  shall  provide  an  estimate  of  these  effects  at  the 
relatively  high  frequency  of  a;  =  30.15  Hz  and  show  that  the  assumption  made  is  indeed  a  reasonable 
one  for  the  frequency  range  considered.  To  arrive  at  an  estimate  of  the  order  of  magnitude,  we  make 
the  simplifying  assumption  that  the  accelaration,  a,  is  constant  and  derive  it  from  the  displacement 
of  the  surface,  a;  =  ^;  thus  The  inertial  force,  Msma^^,  has  to  be  compared  with  the 

actuation  force,  a{t)bd.  We  then  compute  the  ratio  Msma(^/cr{t)bd  =  at  the 

end  of  the  first  half  of  the  first  cycle  (during  the  M  A  transformation)  where  the  actuation  force 
attains  its  maximum.  The  ratio  turns  out  to  be  0.87  %  for  the  case  shown  in  Fig.8;  neglecting 
inertial  effects  thus  appears  to  be  reaisonable, 

7.  CONCLUSIONS 

As  a  first  step  towards  the  design  of  a  high  frequency,  high  force,  large  strain  SMA  actuator, 
we  have  modeled  in  this  work  a  thermoelectrically  cooled  thin  SMA  strip  linear  actuator.  The 
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SMA  undergoes  cyclic  phatse  transition  between  the  martensitic  and  austenitic  phases  by  alternate 
heating/cooling,  achieved  with  the  thermoelectric  Peltier  effect  of  a  pair  of  P/N  semiconductors. 
The  effect  of  variable  actuating  load  and  constant  load  applied  as  boundary  conditions  on  the  SMA 
actuator  have  been  considered.  The  thermomechanical  boundary  value  problem  involves  strongly 
coupled  thermal  and  mechanical  fields.  The  evolution  equations  for  the  field  variables  are  integrated 
using  the  fourth-order  Runge  Kutta  method  and  the  coupling  between  the  fields  is  accounted  for 
by  implementing  an  interative  scheme.  The  primary  parameters  of  interest  in  this  work  are  the 
frequency  response  and  evolution  of  the  variable  load.  The  performance  of  the  actuator  is  compared 
with  various  commercially  available  actuators  based  on  energy  conversion  efficiencies  and  energy 
output  per  unit  volume  of  active  material.  With  reference  to  Fig.8,  we  conclude  that  thin  layers  of 
SMA  about  6  ^  thick,  with  a  low  transformation  hysteresis  and  undergoing  partial  transformation, 
are  potentially  capable  of  delivering  frequencies  of  the  order  of  30  Hz  at  peak  stresses  of  about  145 
Mpa. 


ACKNOWLEDGEMENTS 

The  authors  acknowledge  the  financial  support  of  the  Office  of  Naval  Research,  contract  No.  N 
00014-94-1-0268,  monitored  by  Dr.  Roshdy  Barsoum. 

REFERENCES 

1.  Allen,  D.H,  Thermomechanical  Coupling  in  Inelastic  Solids,  in  Review  of  Thermal  Stresses^ 

Edited  by  Richard  B.  Hetnarski,  Applied  Mechanics  Reviews,  44,  Nos. 8-9,  361  (1991). 

2.  Armero,  F.  and  Simo,  J.C.,  A  New  Unconditionally  Stable  Fractional  Step  Method  for  Nonlinear 

Coupled  Thermomechananical  Problems,  Sudam  Report  No. 91-5,  Stanford  University  (1991). 

3.  Bhattacharyya,  A.,  Lagoudas,  D.C.,  Wang,  A.  and  Kinra,  V.K.,  “On  the  Role  of  Thermoelectric 

Heat  Transfer  in  the  Design  of  SMA  Actuators:  Theoretical  Modeling  and  Experiment”, 5mar^ 
Mater.  Struct.^  4,  252  (1995). 

4.  Bhattacharyya,  A.  and  Lagoudas,  D.C.,  A  Stochastic  Thermodynamic  Model  for  the  Gradual 

Thermal  Transformation  of  SMA  Polycrystals,  Journal  of  Smart  Materials  and  Structures 
(submitted)  (1995). 

5.  Bo,  Z.  and  Lagoudas,  D.C.,  “Deformations  and  Thermal  Response  of  Active  Flexible  Rods 

with  Embedded  SMA  Actuators”,  Smart  Structures  and  Materials  1994:  Smart  Structures 
and  Intelligent  Systems,  Proceedings  of  the  SPIE  1994  North  American  Conference  on  Smart 
Structures,  Orlando,  FA,  2190  (1994). 

6.  Boley,  B.A.  and  Weiner,  J.H., “Theory  of  Thermal  Stresses”,  Krieger  Pub.  Co.,  Florida  (1960). 


28 


7.  Boyd,  J.  and  Lagoudas,  D.C.,  “A  Thermodynamic  Constitutive  Model  for  Shape  Memory 

Materials.  Part  I;  The  Monolithic  Shape  Memory  Alloy  and  Part  II:  The  SMA  Composite 
Material”,  Int.  J.  Plasticity{m -piess)  (1995). 

8.  Brinson,  L.C.,  Bekker  and  Hwang,S.,  “Temperature  Induced  Deformation  in  Shape  Memory 

Alloys”,  in  Active  Materials  and  Smart  Structures,  Edited  by:  Anderson,  G.L.  and  Lagoudas, 
D.C.,  Annual  Technical  Meeting  of  the  Society  of  Engineering  Science,  2427,  234  (1995). 

9.  Buckingham,  E.,  “On  Physically  Similar  System  Illustrations  of  the  Use  of  D3mamical  Equa¬ 

tions”,  Phys.  Rev.,  4,  376  (1914). 

10.  Delaey,  L.,  Krishnan,  R.V.,  Tas,H.  and  Warlimont,H.,  “Thermoelasticity,  Pseudoelcisticity 
and  the  Memory  Effects  associated  with  Martensitic  Transformations”  Parts  I-III  Journal  of 
Materials  Science,  1521  (1974). 

11.  Domenciali,  C.A.,  “Irreversible  Thermodynamics  of  Thermoelectricity”,  Reviews  of  Modem 
Physics,  26,  No.2  ,  237  (1954). 

12.  Giurgiutiu,  V.,  Chaudhry,  Z.  and  Rogers,  C., “Energy-Based  Comparison  of  Solid-State  Ac¬ 
tuators”  ,  Technical  Report,  Center  for  Intelligent  Material  Systems  and  Structures,  VPI-SU, 
Blacksburg,  VA  (1995). 

13.  Harman,  T.C.,  and  Honig,  J.M.,  “Thermoelectric  and  Thermomagnetic  Effects  and  Applica¬ 
tions”,  Mc-Graw  Hill  Book  Co.,  NY  (1967). 

14.  Jackson,  C.M.,  Wagner,  H.J.  and  Wasilewski,  R.J.,  “A  report  on  55-Nitinol  —  The  Alloy  with 
a  Memory:  Its  Physical  Metallurgy,  Properties  and  Applications”,  Technology  Untilization 
Office,  NASA,  Washington,  D.C  (1972). 

15.  Lagoudas,  D.C.,  Bo,  Z.  and  Bhattacharyya,  A.,  “A  Thermodynamic  Constitutive  Model  for 
Gradual  Phase  Transformation  of  SMA  Materials” ,  Proceedings  of  the  SPIE  Conference,  San 
Deigo,  CA,  1996(accepted)  (1996). 

16.  Lagoudas,  D.C.  and  Ding,  Z.,  “Modeling  of  Thermoelectric  Heat  Transfer  in  Shape  Memory 
Alloy  Actuators:  Transient  and  Multiple  Cycle  Solutions” ,”,/nt.  J.  Engng  Sci.,  33,  No. 15,  “ 
2345  (1995). 

17.  Lagoudas,  D.C.  and  Kinra,  V.K.,  “Design  of  High  Frequency  SMA  Actuators,  Disclosure  of 
Invention  Tamus  803,  Texas  A&M  University,  College  Station,  TX  77843  (1993). 

18.  Liang,  C.  and  Rogers,  C.A.,  “One-Dimensional  Thermomechanical  Constitutive  Relations  for 
Shape  Memory  Materials”,  Journal  of  Intelligent  Materials  Systems  and  Structures,  1,  207 
(1990). 


29 


19.  Liang,  C.  and  Rogers,  C.A.,  “The  Multi-dimensional  Constitutive  Relations  of  Shape  Memory 
Alloys” ,  Proceedings,  AIAA  32nd  Structures,  Structural  Dynamics  and  Materials  Conference, 
Baltimore,  MD  (1991). 

20.  Mayergoyz,  I.D.,  “Mathematical  Models  of  Hysteresis” ,  Springer  Verlag,  1983. 

21.  Melcor  Educational  Package,  Melcor,  Trenton,  NJ  (1992). 

22.  Nag,  P.K.,” Engineering  Thermodynamics”,  Tata  McGraw-Hill,  New  Delhi,  India  (1981). 

23.  Numerical  Recipes:  The  Art  of  Scientific  Computing,  Cambridge  University  Press,  London 
(1986). 

24.  Oberaigner,E.R.,  Tanaka,K.  and  Fischer,F.D.,  “Investigation  of  the  Damping  Behavior  of  a 
Vibrating  Shape  Memory  Alloy  Rod  using  a  Micromechanical  Model”,  in  Mathematics  and 
Control  in  Smart  Structures,  Edited  by  V.V.  Varadan,  Proceedings  of  SPIE,  349  (1995). 

25.  Patoor,  E.,Eberhardt,  A.  and  Berveiller,  M.,  “Potentiel  Pseudoelastique  et  Plasticite  de 
Transformation  Martenistique  dans  les  Mono-et  Polycristaux  Metalliques” ,  Acta  Metallurgica, 
35,  No.ll,  2779  (1987). 

26.  Pollock,  D.D.,  “Thermocouples:  Theory  and  Properties”,  CRC  Press,  Boson,  MA  (1991). 

27.  Sun,  Q.P.  and  Hwang,  K.C.,  “Micromechanics  Modeling  for  the  Constitutive  Behavior  of 
Polycrystalline  Shape  Memory  Alloys-I.  Derivation  of  General  Relations”,  J.  Mech.  Phys. 
Solids,  41,  1  (1993a). 

28.  Sun,  Q.P.  and  Hwang,  K.C.,  “Micromechanics  Modeling  for  the  Constitutive  Behavior  of 
Polycrystalline  Shape  Memory  Alloys-I.  Study  of  the  Individual  Phenomena” ,  J.  Mech.  Phys. 
Solids,  41,  19  (1993b). 

29.  Takagi,  T.,  “A  Concept  of  Intelligent  Materials”,  Journal  of  Intell.  Mater.  Sys.  and  Struct., 
1,  149  (1990). 

30.  Taneika,  K.,  “A  Thermomechanical  Sketch  of  Shape  Memory  Effect:  One-Dimensional  Tensile 
Behavior”,  Res  Mechanica,  18,  251  (1986). 

31.  Thrasher,  M.A.,  Shahin,  A.R.,  Meckl,  P.H.,  and  Jones,  J.D.,  “Thermal  Cycling  of  Shape  Mem¬ 
ory  Alloy  Wires  using  Semiconductor  Heat  Pump  Modules”,  Proceedings  of  the  1st  European 
Conference  on  Smart  Structures  and  Materials,  Glasgow,  England  (1992). 

32.  Wada,  B.K.,  Fanson,  J.L.  and  Crawley,  E.F.,  “Adaptive  Structures”,  Journal  of  Intell.  Mater. 

Sys.  and  Struct.,  1,  157  (1990). 


30 


33.  Waymaxi,  C.M.,  “Phase  Transformations,  NondifFusive”,in  Physical  Metallurgy,  Eds:  Cahn, 
R.W.  and  Haasen,  P.  ”,  North-Holland  Physics  Publishing,  NY,  1301  (1983). 

APPENDIX 


A.Determination  of  material  parameters  for  the  thermodynamic  theory 


The  method  of  determination  of  the  parameters  psma^s  ,  psma^u{i),  first  used  in  Eq.(2.12), 
and  Y  in  Eq.(2.18)  from  experiments  are  given. 

Determination  of  psma^s 


The  parameter  psma^s  can  be  found  by  doing  an  uniaxial  tensile  test  at  the  onset  of  transfor¬ 
mation  for  which  ^  =  0.  Prom  Eq.(2.17),  we  have 


_ _ Psma^^ _ 

dT  “  "  {Ell  -  ) AT  -f  H  ’ 


(A.l) 


where  the  total  uniaxial  inelastic  strain,  H,  can  be  determined  from  an  uniaxial  tensile  test  of  the 
polycrystal.  Assuming  that  {E]^  -  -f  -  q;^)AT  <<  H,  Eq.(A.l)  becomes 


d<7  _  PsmgAs 

dT~  H  ■ 

In  the  o  —  T  plane,  a  locii  of  points,  usually  a  straight  line(McCormick  and  Liu,1994),  can  be 
identified  which  represents  the  stresses  at  different  temperatures  for  which  the  material  starts  to 
transform.  Denoting  the  slope  ^  =  D,  we  have 


PsmgAs  —  HD  .  (A. 2) 

While  several  mechanical  tests  are  needed  to  establish  the  line  and  hence  determine  D,  it  is  easier 
to  establish  its  value  from  a  purely  thermal  transformation;  this  is  addressed  shortly. 

Determination  of  psma^uiO 


Psma^u{^)  can  be  determined  from  a  purely  thermal  transformation,  with  <7  =  0  and  <t  =  0. 
The  Eq.(2.11)  reduces  to 


PsmgAs 


Psma 


f  . 


(A.3) 
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where  during  a  purely  thermal  transformation,  we  define 


The  parameter  is 


Psma^^ 


Psma 


(A.4) 


Op  =  —(1  —  ^o)  ^  >  0  and  Op  =  -(^0  C  <  0  , 


(A.5) 


where  ^  =  |o  at  the  onset  of  transformation.  The  function  p{T)  is  taken  to  be  a  normal  distibution, 
given  as(Bhattacharyya  and  Lagoudas,1995) 


p{T)  =  —exp 

bp 


l(T-0MTi  +  Th)\ 

2[  sp  J 


Ti<T<Th, 


.  (A.6) 


where  Sp  is  the  standard  deviation  and 


rn 

=  /  exp 
Jti 


I  T-0.5{Ti+n) 

’2  1  s„ 


dT 


The  function  defined  by  Eq.(A.6)  satisfies  the  following  condition 

r''piT)dT=i. 

Jt, 

The  temperatures  T/  amd  Th  are 

T;  =  MJ  and  Th  =  M^  for  ^  >  0, 

and 

Ti  =  A®  and  for  ^  <  0  , 

where  Mg  and  Mj  are  the  stress-free  martensitic  start  and  finish  temperatures,  and  axe  the 
austenitic  start  and  finish  temperatures.  Using  the  relation  ^  =  app{T)t{ivom  Eqs.(A.3)-(A.4)), 
conceptually  it  is  possible  to  define  T  =  f{^).  Since  for  a  purely  thermal  transformation,  using 
Eq.(A.2),  Eq.(2.17)  reduces  to  H  =  PsmaAsT  -  PsmaAu{0  =  -HDT  —  Psma  Aw  (^),  we  can  write 


Psma^uiO  = -HD f iO +  Y  ,  for  i<0, 


where  Eq.(2.18)  has  been  used.  D  and  Y  are  determined  next.  An  analytical  expression  for  f{^)  is 
not  possible  to  derive  when  p{T)  is  a  normal  distribution;  a  numerical  procedure  is  then  required 
to  determine  /(O  at  any  arbitrary 
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Determination  of  D  and  Y 


For  a  purely  thermal  transformation(a*  =  0  and  &  =  0),  Eqs.(2.15),  (2.17a),  (A. 3)  and  (A.4), 


give 

— V.g  +  r  —  CT  =  ap{psma^sT  —  n)p(T)T  . 
The  total  latent  heat  evolved  is  defined  as 

Hi=  J {-V.q  +  r  -  Ct)dt , 


where 


Hi  ——  OpHa  for  ^<0,Ha>0. 


(A.7) 


(A.8) 


(A.9) 


The  parameter  Ha  is  the  magnitude  of  the  latent  heat  measured  in  a  Differential  Scanning  Calori- 
metric(DSC)  experiment  during  &  A  M  transformation  of  a  SMA  polycrystal  in  an  initially 
austenitic  state(or  vice-versa);  we  have  assumed  for  simplicity  that  the  latent  heat  involved  is 
identical  for  both  transformations.  Eqs.(A.7)-(A.9)  yield 

..  .Ha 


fM°  rA°  1 

1 

i-Mj  fA°f 

/  Tp{T)dT  -  /  Tp{T)dT 

and  Y  =  -HD 

/  Tp{T)dT  -t-  /  Tp{T)dT 

Jm°  Ja° 

2 

_Jm°  Ja% 

(A.IO) 


B.  The  parameters  G{t),  U2{t),  ui  and  F{t) 

and  the  correcponding  nondimensional  parameters 


The  parameters  needed  in  Eq.(5.4)  are 


/i(4,a,r)  =  ^c(0  ^  =  0, 


f 


dAuiO 


HO+ 


dAuio 
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^  (I + s)  /  + 1  /' ■"  f + 


hdd^To  (I  1\  pdd^  ^ 


(A.11) 


whereas  those  in  Eq.(5.23)  are 


G(t)  =  G{t)  ,  m,d,f)  =  /x(e,a,T)  ,  P2(i)  =  . 

p,  =  ,  F(t)  =  ^F{t)  ,  ^(i)  =  Jf(t)  . 


(A.12) 
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Figure  1.  Schematic  diagram  of  a  thin  layer  extensioi 
semiconductor  thermoelectric  elements. 


temperature  and  stress  profiles. 


6(0)  =  3.33  x10- 


Nondimensional  mechanical  work  output,  (b)  energy  efficiencies,  and  (c)  the  corresponding  camot  efficiencies. 


Nondimensional  frequency  of  actuation  at  two  current  levels  for  a  range  of  nondimensional  SMA  layer  thicknesses. 


1.23 


nondimensional  frequency  of  actuation  as  a  function  of  the  hysteresis  in  the  transformation  temperatures. 


=  6  Amps/mm 


Figure  8.  (a)  The  dimensional  temperature  vs.  time  and  (b)  dimensional  stress  vs.  time. 


